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Time-independent searches for astrophysical neutrino sources with the combined data of 40 and 
59 strings of IceCube. 

The IceCube Collaboration 1 

1 See special section in this proceedings 



Abstract: We present the results of time-independent searches for astrophysical neutrino sources performed over the 
whole sky using data collected between April 2008 and May 2010 with the 40-string and 59-string configurations of the 
IceCube Neutrino Observatory. Muon tracks arriving in the detector from neutrino interactions are reconstructed using the 
time and charge information detected by an array of photomultiplier tubes (PMTs). In the northern sky, the data sample 
consists of 14,121 events collected with 40 strings and 43,339 with 59 strings, mostly muons induced by atmospheric 
neutrinos. In this sky region the search is sensitive to point sources of neutrinos with E~ 2 spectra mainly in the TeV-PeV 
energy range. In the opposite hemisphere, a much larger background of high-energy atmospheric muons dominate the 
data set. A zenith dependent energy cut is used to reduce the number of background events. This weakens the sensitivity 
for point sources with E~ 2 spectra with respect to the upgoing region. The downgoing region is more sensitive to harder- 
spectrum sources for which the bulk of events can be detected between PeV-EeV energies. An unbinned maximum 
likelihood ratio test is used to search for astrophysical signals. For the first time it was adapted to combine data from 
different detector configurations. The combined sensitivity is about a factor ~ 2.5 better than the previous 1-year limit of 
the 40-string configuration alone. A dedicated search based on a catalog of sources is also presented. 

Corresponding authors: J. A. Aguilar 2 (aguilar@icecube.wisc.edu), M. Baker 2 , J. Dumm 2 , T. Montaruli 2 , 
N. Kurahashi 2 

2 Dept. of Physics, University of Wisconsin, Madison, WI 53706, USA 
Keywords: neutrino astronomy; point sources; likelihood method 



1 Introduction 

The IceCube Neutrino Observatory is a neutrino telescope 
installed in the deep ice at the geographic South Pole. The 
final configuration comprises 5,160 photomultiplier tubes 
(PMTs) [1] along 86 strings instrumented between 1.5-2.5 
km in the ice sheet. Its design is optimized for the detec- 
tion of high energy astrophysical neutrinos with energies 
above <~ 100 GeV. The observation of cosmic neutrinos 
will be a direct proof of hadronic particle acceleration and 
will reveal the origins of cosmic rays (CR) and the possi- 
ble connection to shock acceleration in Supernova Rem- 
nants (SNR), Active Galactic Nuclei (AGN) or Gamma 
Ray Bursts (GRBs). The IceCube detector uses the Antarc- 
tic ice as the detection volume where muon neutrino inter- 
actions produce muons that induce Cherenkov light. The 
light propagates through the transparent medium and can 
be collected by PMTs housed inside Digital Optical Mod- 
ules (DOMs). The DOMs are spherical, pressure resis- 
tant glass vessels each containing a 25 cm diameter Hama- 
matsu photomultiplier and its associated electronics. De- 
tector construction finished during the austral summer of 
2010-11. 



During its construction, the IceCube telescope ran in var- 
ious configurations. From April 2008 to May 2009, 40 
strings of IceCube were operational and collecting scien- 
tific data. The data analyzed for that period has provided 
until now the best sensitivity to high energy neutrino point 
sources. In this article we describe the analysis of the com- 
bined data of the previous 40 strings of IceCube with the 
59-string configuration data from May 2009 to May 2010 
(see Fig. 1). With the combined information from both data 
samples we are able for the first time to probe beyond the 
initial estimates of the 1-year sensitivity of the completed 
IceCube detector [2]. 

2 Method 

The method used for this analysis is an unbinned maxi- 
mum likelihood ratio test [4]. This test allows to calcu- 
late the significance of an excess of neutrinos above the 
background for a given direction. The method uses both 
the reconstructed direction of the events as well as the re- 
constructed visible muon energy, to discriminate between 
signal and background [3]. This method improves the sen- 
sitivity to astrophysical sources over directional clustering 
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3 Data selection and Detector Performance 
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Figure 1: Detector layout. The empty circles represents 
the string positions that corresponds to the geometry of the 
whole IceCube detector. The 40-string configuration is rep- 
resented with the small dots and the 59-string configuration 
with filled circles. 



alone by leveraging the event energies in order to separate 
hard-spectrum signals from the softer spectra of the atmo- 
spheric neutrino or muon backgrounds. For each tested di- 
rection in the sky, the number of signal events over back- 
ground, n s , and the index of a power law, 7, for the spec- 
trum of the signal events are determined that maximize the 
likelihood function. The likelihood ratio between the best- 
fit hypothesis and the null hypothesis (n s = 0) forms the 
test statistic. The significance of the result is evaluated by 
performing the analysis on scrambled data sets, randomiz- 
ing the events in right ascension but keeping all other event 
properties fixed. Uniform exposure in right ascension is 
ensured by the daily rotation of the detector with respect 
to the sky. Events that are close to the polar regions of the 
sky (declination < —85° or > 85°) are excluded from the 
analysis, since the scrambling in right ascension does not 
work in these regions. 

Two point-source searches are performed. The first is an 
all-sky search where the maximum likelihood ratio is eval- 
uated for each direction in the sky on a grid of 0.1° x 0.1°, 
much finer than the angular resolution. The significance 
of any point on the grid is determined by the fraction of 
scrambled data sets containing at least one grid point with 
a likelihood ratio higher than the one observed in the data. 
This fraction is the post-trial p-value for the all-sky search. 
Because the all-sky search includes a large number of ef- 
fective trials, the second search is restricted to the direc- 
tions of a priori selected sources of interest. The post-trial 
p-value for this search is again calculated by performing 
the same analysis on scrambled data sets. 



From April 2008 to May 2009, 40 strings of the IceCube 
detector were operational. The duty cycle at analysis level 
for that period was ~ 90% after selecting good runs based 
on the detector stability and the total livetime was 375.5 
days. The event selection for the 40-string configuration 
data sample is described in detail in Ref. [3]. In the analy- 
sis of the data from the 40 strings of IceCube no significant 
excess over fluctuations of the background was found, and 
upper limits have been published (see Ref. [3]). Here we 
report on a combined analysis using the data correspond- 
ing to the 40 strings period plus the 59-string configuration 
data taken from May 2009 to May 2010. The livetime cor- 
responding to the 59-string configuration is 348 days with 
a similar ~ 90% duty cycle at final analysis level as in the 
previous year. 

The trigger rate of the 59-string configuration is of the or- 
der of 1.5 kHz for events based on a simple multiplicity 
trigger requiring 8 triggered DOMs. This trigger rate is 
strongly dominated by the muon background produced in 
cosmic rays interactions in the atmosphere. A first level 
of background rejection of poorly reconstructed up-going 
events and a selection of high energy muons for the south- 
ern sky is done on-site at the South Pole (Level 1 filter). 
The data sent through the satellite to the North undergo 
further processing that includes a broader range of more 
CPU consuming likelihood-based reconstructions at the so- 
called Level 2 and Level 3 filters. This offline processing 
also provides useful parameters for background rejection 
and measurements of the energy and of the angular uncer- 
tainty. The data rate at Level 3 is of the order of 3 Hz and 
still dominated by atmospheric muons. Because the north- 
ern sky and the southern sky present very different chal- 
lenges, two separate techniques for background rejection 
are used for each hemisphere. 

In the northern sky the 59-string configuration event se- 
lection was performed using a multi-variate classification 
algorithm. Boosted Decision Trees (BDTs) were used in 
the final analysis step to classify events as signal-like or 
background-like. Twelve event observables, split in two 
sets of 8 and 4 respectively, were selected by choosing 
variables with low correlation for a background dominated 
dataset (correlation coefficient |c| < 0.5), but high discrim- 
inating power between signal and background. Training 
was done using a subsample of the data as the background 
and simulated neutrino events as signal. We trained two 
sets of BDTs, one with a neutrino spectrum of E~ 2 , and 
one with a neutrino spectrum of E~ 2J in order to account 
for softer neutrino spectra and possible TeV cut-offs in the 
expected neutrino emission. The usage of data as the back- 
ground sample for training is an important aspect since it 
makes the analysis independent of the systematic uncer- 
tainties of the simulation of the muon background. The 
combination of the two BDT scores for both the softer neu- 
trino spectrum and the standard E~ 2 is used as a selection 
criterion and the cut that optimizes the discovery potential 
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for neutrino point sources over a wide energy range is cho- 
sen. 

The southern sky is filtered by using energy estimators 
to separate the large amount of down-going atmospheric 
muons from a hypothetical neutrino signal with a harder 
spectrum. For vertically down-going events with zenith an- 
gles between and 50 degrees, we take advantage of the 
IceTop detector in order to reject atmospheric muons orig- 
inating from a shower that produces a signal in at least two 
of the PMTs of the IceTop detector. This IceTop veto al- 
lows us to reject background with ~ 99% efficiency in the 
very vertical zenith angles without losing signal neutrino 
efficiency (< 1 %). 
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Figure 2: Solid angle averaged neutrino effective area in the 
northern sky for the 59-string IceCube configuration (dot- 
ted line) and the 40-string configuration (solid line) for an 
equal ratio of and as a function of the true neutrino 
energy. 
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Figure 3: Nuetrino angular resolution as a function the 
true neutrino energy for the 59-string IceCube configura- 
tion and the 40-string configuration. 



Figure 4: Expected sensitivity (solid line) for 90 % C.L. 
using the classical (frequentist) construction of upper lim- 
its, and discovery potential defined as the minimum flux 
required to have a 50% probability to claim a discovery 
of a point-source with a E~ 2 neutrino spectrum with con- 
fidence level equivalent to 5a, (dotted line), for the com- 
bined analysis using the 40-string and 59-string configura- 
tion data. Both lines are shown as a function of the decli- 
nation. 

come from the southern sky and the rest from the up-going 
region. The estimated atmospheric muon contamination in 
the northern sky is ~5 %. The solid angle averaged neu- 
trino effective area for both detector configurations in the 
northern sky is shown in Fig. 2 as a function of the true 
neutrino energy. The overall increase in neutrino effective 
area in the up-going region of the 59-string configuration 
with respect to the previous IceCube configuration of 40 
strings is a factor of ~ 1.3 for energies > 1 TeV and up 
to a factor of 2 at lower energies due to the event selection 
based on BDTs trained with softer neutrino spectra. Fig. 3 
shows the angular resolution defined as the median of the 
point spread function (PSF) as a function of the true neu- 
trino energy. The PSF is defined as the angle between the 
reconstructed muon track and the true neutrino direction. 
The BDT used in the 59-string configuration allows more 
low energy signal events to pass the event selection with 
worse angular resolution, that explains why the median dis- 
tribution is worse compared to the 40-string configuration 
at energies below 10 TeV. 

The expected sensitivity for the 2 years (375 + 348 days) 
of combined data and the discovery potential is shown in 
Fig. 4 as a function of declination for a E~ 2 neutrino spec- 
trum. The overall improvement with respect to the 40- 
string configuration sensitivity is about a factor of ~ 2.5 
making it comparable to the projected 1-year sensitivity of 
the completed IceCube detector. 



The final data sample for the 59-string configuration has 
a total number of 107,569 events, among them almost 2/3 
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Figure 5: Significance skymap in equatorial coordinates (J2000) of the all-sky point source scan. The dotted line indicates 
the galactic plane. 



4 Results 

The results of the all-sky scan are shown in the preliminary 
pre-trial significance map in Fig. 5. The most significant 
deviation from background is located at 75.45° r.a. and 
— 18.15° dec. The best-fit parameters for this location are 
fis = 18 and 7 = 3.9. The pre-trial p-value is 2.23 x 10 -5 
which corresponds to a post-trial p-value of 67% calcu- 
lated as the fraction of scrambled sky maps with at least 
one source with an equal or higher significance. The most 
significant source from the a priori source list is PKS 1454- 
354 with a pre-trial estimated p-value of 14%. The equiv- 
alent post-trial p- value was calculated as well using scram- 
bled sky maps and correspond to a value of ~ 95%. The 
90% CL upper limits for both searches will be provided as 
the systematic uncertainties are evaluated and incorporated 
in the Feldman & Cousins confidence belt construction. 



of this analysis is already beyond the initial estimates of 
the expected sensitivity of 1 year of the IceCube detec- 
tor. Two searches were performed; an all-sky scan and a 
search on specific locations based on a priori list of candi- 
date sources. In both cases no significant excess was found. 
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5 Conclusions 

Between April 2008 and May 2009 the IceCube detector 
recorded 375 days of data with 40 instrumented strings. 
The analysis included 36,900 events in the whole sky where 
no evidence for a signal was found. The 40-string configu- 
ration analysis provided the best flux upper limits on point 
sources of astrophysical neutrinos up to now. Here we pre- 
sented an analysis of a combined data sample in 40-string 
and 59-string configurations of IceCube. The sensitivity 
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Searches for Time- Variable Neutrino Point Sources with the IceCube Observatory 

The IceCube Collaboration 1 
1 See special section of these proceedings. 



Abstract: We present searches for time-dependent emissions of neutrinos in the entire sky using the data collected 
between April 2008 and May 2010 with 40 and 59 strings of IceCube. An all-sky search is performed searching for 
any clustering of events in space and time. In the northern sky the sample is mainly atmospheric neutrinos, while 
in the southern sky the sample is dominated by atmospheric muons. In order to reduce the penalty of trials we also 
perform a search based on flares of AGNs observed by other experiments, using lightcurve information from bands where 
comprehensive coverage is available. Results from the 40-string detector are presented in this paper, while those from the 
59 string detector will be presented at the conference. 
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1 Introduction 

The IceCube Neutrino Observatory is a kilometer-scale de- 
tector located at the geographic South Pole. Beneath the 
glacial surface, IceCube is composed of 5160 optical mod- 
ules (DOMs) deployed on 86 vertical strings between 1450 
and 2450 m to detect and reconstruct high energy neutrino- 
induced charged leptons. The main science goal of the Ice- 
Cube experiment is the detection of astrophysical neutrino 
sources, which will help identify the origins of the highest 
energy cosmic rays. 

Muons passing through the detector emit Cerenkov light 
allowing reconstruction with median angular resolution of 
less than 1° for > 10 TeV energy muons in the 40 and 59- 
string configurations. 

Time-dependent analyses aim to reduce the background 
of atmospheric neutrinos and muons by searching over 
smaller time scales around a period of interest. The 
searches discussed in this paper are about a factor of four 
more powerful than time-integrated searches for flares with 
duration < 100 s. In this paper we describe the addi- 
tion of a time dependent term to the standard searches for 
steady emission of neutrinos presented in [1] [2]. We ap- 
ply this term in searches for neutrino emission from an all- 
sky generic time-dependent search and from a catalogue 
of sources with photon flares occurring when IceCube was 
taking data in its 40 and 59-string configurations. 



2 Time Dependent Point Source Searches 

An unbinned maximum likelihood ratio method which 
models the data as a mixture of signal and background has 
been used for the search for point sources of neutrinos in 
IceCube [3]. We use the angular and energy distribution 
of events as information to characterize the signal with re- 
spect to the background. Astrophysical sources of neutri- 
nos will cluster near the object and are expected to have a 
power-law energy spectrum harder than that measured of 
the atmospheric backgrounds [4]. An energy estimator is 
used based on the photon density along each reconstructed 
muon track. The atmospheric background has a roughly 
constant rate in time, but sources such as Active Galactic 
Nuclei exhibit significant variability in photon flux states, 
allowing for tests aimed at additional background rejection. 
The analysis returns a best-fit number of signal events and 
spectral index, as well as other free parameters from the 
time-dependent terms. 

The IceCube 40-string data at analysis level consists of 
36,900 selected events, 14,121 are upward-going neutrino 
candidate events and 22,779 are downward-going, mainly 
PeV energy muons from atmospheric air showers [1]. The 
59-string data at analysis level consists of 107,569 selected 
events, 43,339 of which are upward-going neutrino candi- 
dates, and 64,230 are in the downward-going region [2]. 
The data used in this analysis were collected over 724 days 
of livetime between April 5, 2008 and May 31, 2010, cor- 
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responding to 92% of all data taken during that period. Se- 
lection cuts for the final sample are based on the quality 
of the reconstruction, such as the angular uncertainty of 
the track reconstruction, and on other variables such as the 
number of DOMs hit by the direct Cerenkov light produced 
by muons. These variables help reject misreconstructed 
events. 

The signal probability distribution function (pdf) is: 

Si = Sf ace (l Si - x s \,Oi)Sr erav {Ei,6 u ls)S? me , 

(1) 

where cr^ is the reconstructed angular uncertainty of the 
event [5], | x*i — x s | the angular separation between the re- 
constructed event and the source, Sl ner9V (E i7 9i,j s ) is the 
energy pdf with the event energy Ei, zenith angle 6i, and 
spectral index 7 S , which is built in 22 zenith bands to ac- 
count for declination dependence of the background. Sf me 
is the time pdf of the event. The background pdf is given 
by: 

Bi = B 8 i pace (9i,4>i) Bi ner9V (Ei,9i) Bf me (U), (2) 

where Bl pace (6 i ,4> i ) is the background event density (a 
function of the azimuth <j>i and zenith 9i of the event), 
Bi nergv (Ei, 9i) is the zenith-dependent energy distribution 
of the background, and Bj me (ti) the inverse of the live- 
time, reflecting the fact that the probability density func- 
tions are normalized to one and the background rate is es- 
sentially fiat in time. The background pdf is determined 
using the data, and the final p-value for each analysis is ob- 
tained by comparing equivalent experiments scrambled in 
time and right ascension to actually observed data. 

3 All-Sky Time Scan 

We perform a scan for any significant excess with respect to 
background over all time scales (from sub-seconds to a full 
year) over the entire sky. Since this analysis finds events 
clustered in time, independent scans are performed using 
the 40-string (April 5, 2008 to May 20, 2009) and 59-string 
(May 20, 2009 to May 3 1, 2010) samples. For flares shorter 
than ^100 days, this provides a better discovery potential 
for the time-dependent hypothesis than a time-integrated 
analysis. In principle short bursts can be discovered at a 5a 
threshold with only two events if they occur near enough 
in time. An advantage of an untriggered search such as 
this is the ability to probe all emission scenarios, including 
neutrino emission without any observed counterpart in the 
electromagnetic spectrum. 

This analysis method was developed and tested using a 
simulation of a generic km 3 neutrino detector in [6], and 
has been adapted for use with a detector with non-uniform 
acceptance in zenith and azimuth and dead time [7]. The 
time-dependent probability density function for this search 
is a Gaussian function, with its mean and width as free pa- 
rameters, returning the most significant flare from a partic- 
ular direction. The method is applied as an all-sky scan 



over a 0.5°x 0.5° grid (smaller than the typical event an- 
gular uncertainty) in right ascension and declination, scan- 
ning for flares of all durations from 20 ^s (the minimum 
time separation between events) to 150 days. The final re- 
sult is the set of best fit parameters from the location with 
the strongest deviation from background. 
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Figure 1: The 50% 5er discovery potential and 90% sen- 
sitivity in terms of the mean number of events for a fixed 
source at +16° declination with the 59-string detector. The 
number of events for the sensitivity and discovery poten- 
tial for the time-independent search are also shown. Flares 
with a a T of less than 40 days, or a FWHM of less than 
roughly 100 days during the 59-string data taking period, 
have a better discovery potential than the steady search. 



3.1 40-String Results 

Using the 40-string data, the location which deviates 
most from background is found at (RA,Dec)=(254.75°, 
+36.25°), and is presented in [7]. Two events are found 
(2.0° apart in space and 22 seconds in time), with a best- 
fit spectrum 7 of 2.15 (with uncertainty of ±0.4), mean of 
the flare T a of MJD 54874.703125 and width a T of 15 sec- 
onds. A clustering of greater significance is seen in 56% of 
scrambled skymaps, which is consistent with a fluctuation 
of the background. The result of the scan with 59-strings 
data will be presented at the conference. 

4 Triggered Search for Flares 

When there is specific timing information about the pho- 
ton activity of an astronomical object, that information can 
be used to reduce the background. For triggered sources, 
the focus is on objects such as blazars, which exhibit vari- 
ability on timescales of hours to weeks. When flares are 
seen with comprehensive coverage, flux measurements are 
made on a regular basis and this continuous lightcurve can 
be used to define the activity at any point in time as low 
to high. This improves the ability to define periods of high 
flux state with a clear beginning and end. The assumption is 
that photon and neutrino emissions experience heightened 
states simultaneously. This analysis utilizes 1-day binned 
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lightcurves from the Fermi LAT. Results for 40-string data 
are presented in [7], while at the conference results includ- 
ing the 59 string period will be shown. 

4.1 Method and Expected Performance 

Sources for this search were selected considering alerts 
during the 40 and 59-string data taking periods for sources 
in outburst > 1.5 x 10~ 6 photons/s/cm 2 . Sources with 
flares in the 40-string period (April 5, 2008 to May 20, 
2009) are listed in table 1 . A Maximum Likelihood Block 
(MLB) algorithm [8] [9] is used to de-noise the lightcurves 
by iterating over the data points to select periods from the 
lightcurves which are consistent with a constant flux, tak- 
ing statistical errors into account. The hypothesis is that the 
neutrino emission follows the lightcurve, but only when the 
photon flux goes above a certain threshold F th . By look- 
ing only at these high states the atmospheric background is 
largely reduced. The value of F t h is used as a free parame- 
ter, finding the value of the threshold which maximizes the 
significance of the data. This method is designed to avoid 
any penalty from making an incorrect a priori choice on 
a flaring threshold, which is larger than the effect of one 
additional degree of freedom. 

F(ti) is defined as the value of the denoised light curve at 
ti and F t h is the flux threshold below which no neutrino 
emission is assumed (i.e. Sf me =0 if F(t t ) < F th ). In the 
case of F(ti) > F t h, the probability of neutrino emission 
is assumed to be proportional to the flux level above that 
threshold: 

otime _ ( F (ti) - F th ) , 

1 ~ Nf ' { ' 

where the normalization factor Nf is the integral of the de- 
noised light curve above the threshold. Allowing a lead 
or lag of up to 50 days was also tested. This resulted in 
a markedly higher number of events for discovery, so we 
constrained the neutrinos to come within ± 1 day of the 
photons. 



Source 


pre-trial 


Threshold (10~ 6 


Duration above 




p-value 


cm -2 s _1 ) 


threshold (days) 


PKS 1510-089 


— 





282 


3C 66A/B 


0.47 


0.675 


57 


3C 454.3 


0.20 


9.47 


2.5 


PKS 1454-354 







282 


3C 279 


0.47 


2.34 


6 


PKS 0454-234 







282 


PKS 1502+106 


0.049 


3.13 


8 


J 123939+044409 







282 



Table 1 : Sources tested with the 40 string data and pre-trial p- 
values for the flare search with continuous lightcurves. In the 
event of an underfluctuation no p-value is calculated. The overlap 
between the Fermi public release data and the 40-string data tak- 
ing period is 282 days, that being the maximum duration of the 
lightcurve above F th . 



most significant cluster of events is separated in time by 22 
s and in space by 2° and has a trial-corrected p-value of 
56%. The most significant observation of a flare from cata- 
logues compiled using Fermi-LAT and IACT alerts during 
the 40-string configuration data taking is PKS 1502+106, 
with a p-value of 29% after trials. The results of time- 
dependent searches including the IceCube 59-string data 
will be presented at the conference. 
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4.2 Results 



The results from all sources tested during the 40-string 
data-taking are listed in table 1 . The most significant source 
from the 40-string data-taking period is PKS 1502+106, 
which has a pre-trial p-value of 5%. With the method, we 
find one high-energy event during the August 2008 flare. 
The post-trial p-value is 29%, which is compatible with 
background fluctuations. Results extending the lightcurves 
and adding the 59-strings data will be presented at the con- 
ference. 



5 Conclusion 

We have analyzed data from the IceCube observatory 
from the season 2008-9 when the detector consisted of 40 
strings. The all-sky scan over all directions finds that the 
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Figure 2: An example of the one-day binned Fermi lightcurve (gray points, with statistical errors) and denoised lightcurve 
(light gray solid line) for the blazar 3C273. The lightcurve begins here when Fermi science operations began, the time 
axis continues until the end of the 59-string data taking period. 
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Abstract: The combined data of IceCube's 40-string and 59-string configurations spanning over 700 days are used to 
search for astrophysical neutrinos originating from sources listed in several catalogs. A stacking method which stacks 
sources and searches for an integrated signal above the estimated background is employed. Very large scale sources such 
as neutrino emission from the galactic plane and the Fermi bubble region are also used as extended source hypotheses. To 
perform these searches, a likelihood method that tests the presence of signal using the shape of such sources, their energy 
spectra and the angular distribution of events is employed. 
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1 Introduction 

A search for astrophysical neutrinos originating in galactic 
and extragalactic sources using a likelihood (LH) method 
extensively described in [1, 2] has been performed re- 
cently on IceCube data and is reported in [3]. This method 
uses energy and directional information that distinguish the 
softer backgrounds of atmospheric muons and neutrinos 
from the harder astrophysical neutrinos. Such neutrinos, 
yet to be observed, could originate in jets and shocks via 
Fermi acceleration. The background is estimated by scram- 
bling real events in their arrival times (or right ascensions) 
in the LH method. In this way, the p-value (the fraction of 
randomized data sets with higher test-statistic values than 
the real data) comes only from data and has no dependency 
on the accuracy of the simulation [1]. While Ref. [3] is 
focused on the LH search for steady emissions from point 
sources from the whole sky and from selected sources of 
interest, the work reported in [4] and [5] is focused on 
extending the likelihood method by utilizing time depen- 
dence of emissions (GRBs and AGN flares). The work 
presented here extends the LH method to the stacking of 
sources belonging to the same source class and to very ex- 
tended sources that cover a large fraction of the sky. In 
this work, as in [1, 3], the median sensitivity and upper 
limits at a 90% confidence level (CL) are calculated ac- 
cording to teh classical (frequentist) construction of upper 
limits (Neyman 1937). The discovery potential is the flux 
required for 50% of trials with simulated signal to yield 
a p-value less than 2.87 x 10~ 7 (i.e. 5a significance if 



expressed as the one-sided tail of a Gaussian distribution) 
unless stated otherwise. 

The data sample and its selection is described in [3]. For 
the first time the LH method has been adapted to combine 
data belonging to different configurations. We combine the 
data sample collected for 375.5 d in the 40-string configu- 
ration of IceCube during the period from 2008 April 5 to 

2009 May 20 [1] and data collected during the 59-string 
configuration for 348.1 d from May 20, 2009 to May 30, 

2010 [3]. The total data sample consists of 43,339 (64,230) 
events with 59 strings in the upgoing (downgoing) hemi- 
sphere and 14,121 (22,779) for 40 strings. Hence the total 
number of events on which this search is performed in the 
whole sky is 144,469. As shown in [3], depending on the 
declination the sensitivity of the 40+59-string sample with 
respect to the 40-string sample published in [1] improves 
by about a factor of 2.5. A factor of about 4 is achieved 
in the downgoing vertical region where use of the IceTop 
detector as a veto for muons leads to a more significant 
improvement [3]. It is to be noted that the search for astro- 
physical sources in IceCube extends to the entire sky but 
the sensitivity is different in the upgoing (Northern sky) 
and in the downgoing (Southern sky) regions. As a matter 
of fact, as explained in [1], the upgoing region is domi- 
nated by atmospheric neutrinos since muons are filtered by 
the Earth, while we select high energy muons in the down- 
going region since we look for a clustering of astrophysical 
events characterized by a harder spectrum than atmospheric 
events. 
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All searches shown in this proceeding will be unblinded by 
the time of the conference, hence here we indicate prelim- 
inary sensitivities. P- values and corresponding fluxes will 
be given at the conference. 

2 Stacking searches 

Stacking multiple sources in neutrino astronomy can en- 
hance discovery potential and further constrain astrophysi- 
cal models for uniform populations of sources. The stack- 
ing method is described in detail in [1] where it is explained 
how the signal and background are integrated over a set of 
sources using the same weight for all sources or a weight- 
ing scheme from models for specific tests. As shown in [1], 
the fractional flux needed for discovery for stacked sources 
compared to single sources at 5a CL is very close to the in- 
verse of number of stacked sources. The stacking searches 
we perform were also performed in [1] on the 40-string 
sample with the exception of one catalog. We perform: 

1. A stacking search for 17 Milagro TeV gamma 
ray sources, 9 of which have been reported as 
high significance-detections and another 8 which are 
lower in significance but also confirmed by Fermi [6, 
7]. The sources include supernova remnants, pulsars, 
pulsar wind nebulae, and one unconfirmed hot-spot; 

2. It was noticed with an a posteriori search [1], that 
some of the 6 supernova remnants (SNR) observed 
by Milagro at energies > ITeV and considered in 
[8] as interesting potential neutrino emitters are the 
most significant in the previous list above. When 
these 6 sources are analyzed as a single sub-group, 
an a posteriori p-value of 0.02 was found with best 
fit parameters n s = number of signal events =15.2 
and 7 = spectral index = 2.9. The true trial factor 
is incalculable from this search with 40 strings since 
this was done after unblinding. Hence we perform a 
search for these 6 sources using only the sample of 
59-strings in order not to bias the final result; 

3. a stacking search for 127 local starburst galaxies [9]; 

4. a stacking search for five nearby galaxy clusters 
(GCs), testing four different models for the CR spa- 
tial distribution [10]. The GCs we considered are 
Virgo, Centaurus, Perseus, Coma and Ophiuchus. 
The parameterization of the models, taken from [10] 
and described in detail in [1], consider scenarios 
where CRs are uniformly distributed within the clus- 
ter shock radius or the virial radius (Models A and B 
respectively), an Isobaric Model in which CRs fol- 
low the distribution of thermal gas, and a Central 
AGN Model in which CRs are accelerated in a cen- 
tral AGN. This last case is treated as a point source 
while other models include the extension of sources 
in the LH. 



5. a stacking search of ultra-high energy cosmic rays 
(UHECRs) assuming the sources of these UHECRs 
are also neutrino sources. We include UHECRs ob- 
served by the Auger Observatory [11] and the HiReS 
[12] collaboration. 

The UHECR search presented here looks for neutrinos 
originating from the same direction as observed UHE- 
CRs. However, for this search, we fit for the 'extension' 
of the sources in maximizing the LH. The 'extension' of 
the sources here represent the possible deflection of the 
UHECR from their sources which depends on the UHE- 
CRs energy and the distance and magnetic field strength 
they propagate through. The observation of a correlation 
between UHECR and neutrinos is only possible in proton 
dominated scenarios since heavy elements would be de- 
flected too much by magnetic fields. This search fits for an 
average deflection of all events stacked. The average can 
still be insightful given the similar energies of the UHE- 
CRs. Furthermore, the width of a Gaussian centered at the 
UHECR directions is used as the fit variable, accounting for 
the possible event-by-event directional variation between 
the UHECRs and neutrinos. This search, like all other LH 
searches performed, incorporates the neutrino event point 
spread function of the IceCube detector as well. After fit- 
ting for the 'extension' of sources, we test a model describ- 
ing the deflections of UHECRs in extragalactic magnetic 
fields [13] and constrain the degenerate parameter space 
of distance and B-field intensity assuming that UHECRs 
are protons. The Auger sample contains 69 events with 
reconstructed energy above 55 EeV [11] with angular res- 
olution better than 0.9° for events that trigger more than 6 
stations, and the HiRes sample contains 13 events [12] in 
the same energy range with angular resolution of 0.8° in 
stereo mode. The energy scale for these events is known at 
the level of 20%. In order to reduce the galactic magnetic 
field contribution, expected to be larger than the intergalac- 
tic one, UHECR events that cross the galactic plane are 
removed so the fitted extension represents the intergalactic 
magnetic deflection alone. 

The sensitivities for these searches are shown in Tab. 1 . 

3 The Galactic plane 

Neutrinos and photons are expected to be produced in the 
galactic plane via interaction of cosmic rays with the inter- 
stellar medium (ISM). Diffused cosmic rays are confined 
for million years in the galaxy and their spectrum is hence 
expected to approach E~ 2 - 7 with a composition that be- 
comes heavier above the knee at a few PeVs. The low den- 
sity of the ISM allows secondary mesons to decay before 
reaching their interaction length, thus preserving the cos- 
mic ray spectrum in their decay products and also fixing 
the ratio of neutrinos to photons. The Milagro collabora- 
tion have reported observation of TeV gamma-rays from 
the galactic plane [19], showing good agreement with the 
GALPROP model. The GALPROP code propagates cos- 
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Catalog 


Data sample 


Fit par. 


Sensitivity 


Disc. Pot. 


Milagro 17 


I+II 


l+n s 


0.28 x E2 


1.05 X.B2 


Milagro 6 


II 


l+n s 


0.66 x FM 


3.20 xFM 


Starbursts 127 


I+II 


l+n s 


0.96 x E2 


3.44 xE2 


GC - Model A 


I+II 


l+n s 


1.39 x FM 


4.85 xFM 


GC - Model B 


I+II 


l+n s 


2.58 x FM 


8.98 xFM 


GC - Isobaric 


I+II 


l+n s 


1.26 x FM 


4.84 xFM 


GC - AGN 


I+II 


l+n s 


0.63 x FM 


2.36 xFM 


UHECR 


I+II 


•y+a+n s 


4.01 x E2 


12.8 xE2* 



Table 1: Median sensitivities of the stacking searches (90% CL) and the discovery potentials (p-value < 2.87 x 10 -7 ) 
given as a sum of flux required from the sources in each catalog. Data sample I (II) represents data taken with the 40-string 
(59-string) configuration. Fit parameters are the spectral index 7, the extension of the sources a, and the number of signal 
events n s . E2 indicate units of E 2 dN/dE [10 -11 TeV cm" 2 s" 1 ], while FM indicates the scaling to the predicted flux 
profile. *For UHECR, p-value < 1.35 x 10~ 3 (3a significance of one-sided Gaussian) is used for discovery potential 
calculation. 




mic rays in the galaxy with assumptions of the distribution 
of cosmic ray sources [15]. The Fermi telescope has prod- 
cued a detailed map of the observed gamma-ray emission 
from the galactic plane region [17]. The GALPROP code 
and a fit for the cosmic ray flux are then used to separate the 
contribution from neutral pion decays, which the neutrino 
emission map should closely follow. We use this Fermi 
model of the relative flux of pion decays as the spacial 
template in the LH analysis for neutrinos from the galac- 
tic plane. Fig. 1 shows the relative neutrino signal strength 
expected in the 59-string IceCube configuration assuming a 
source emission pattern of the Fermi model. The strongest 
emission is expected near the galactic center and extends 
along the plane. However, because the center belongs to 
the Southern (downgoing) hemisphere, the signal strength 
is expected to be highest at the closest region on the plane 
to the galactic center that belongs to the Northern (upgoing) 
hemisphere, due to IceCube's background and event selec- 
tion as described in the Introduction section. The LH anal- 
ysis will combine data from both the 40 and 59-string con- 
figurations. One noteworthy issue is that while the Fermi 
model used here provides a relative flux expected from dif- 
ferent regions of the sky, it does not provide the absolute 
scale of the predicted flux. Therefore, a detailed calcula- 
tion of the neutrino flux must be made, or older flux predic- 
tions [16] must be used to translate event counts in IceCube 
into fluxes. The flux calculations of [16] predict 21 neu- 
trino events in the combined 40-string and 59-string data. 

4 The Fermi Bubbles 

Recently evidence has emerged of enormous features in 
the 7-ray sky observed by the Fermi-LAT instrument: bi- 
lateral "bubbles" of emission centered on the core of the 
Galaxy and extending to around 10 kpc from the Galactic 
plane. These structures are coincident with a nonthermal 
microwave "haze"' and an extended region of X-ray emis- 
sion [18]. The bubbles 7-ray emission is characterized by 
a hard and relatively uniform spectrum, relatively uniform 
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Figure 1: Relative signal strength from different parts of 
the sky expected from the galactic plane to be seen with 
the 59-string IceCube detector, plotted in galactic coordi- 
nates (preliminary). The center belongs in the Southern 
(downgoing) hemisphere, so the signal strength is highest 
at the closest region on the plane to the galactic center that 
belongs to the Northern (upgoing) hemisphere as expected. 

intensity, and an overall luminosity of <~ 4 x 10 37 erg/s. 
The 7-ray luminosity between 1-100 GeV is measured to 
be around an order of magnitude larger than the microwave 
luminosity and more than an order of magnitude lower than 
the X-ray luminosity. In [20] the bubbles are explained as 
due to a population of relic cosmic ray protons and heav- 
ier ions in the energy range of 10 — 1000 GeV trapped for 
timescales of about 10 10 yrs undergoing pp collisions on 
the bubbles low density plasma that produce secondaries 
which in turn produce 7-rays and neutrinos. Accounting 
for ionization losses by sub-relativistic protons and adia- 
batic energy losses at all energies, bubble protons lose a 
total of about 10 39 erg/s in steady state, precisely account- 
ing for the CR power injected at the galactic center. This 
is compatible with observations in 7-rays around the galac- 
tic center by Fermi and at higher energies by HESS [20]. 
While in the galactic plane diffusive confinement of the 
CRs leads to a steepening of the steady-state spectrum to 
oc E~ 2 - 7 , in the bubbles there is no energy-dependent con- 
finement effect. So, given the almost energy-independent 
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Figure 2: Relative signal strength expected from Fermi 
bubbles with the 59-string IceCube detector, plotted in 
galactic coordinates (preliminary). Most of the signal is 
expected from the corner of the bubble that is in the North- 
ern (upgoing) hemisphere when transformed to equatorial 
coordinates. 

pp loss time, the observed spectrum of the CRs is as at 
the injection at their acceleration sites (evidently oc E~ 21 ) 
mirrored by the bubble 7— rays. If the bubbles are a reser- 
voir of CRs it is possible to think that the primary spectra 
reach energies above the knee and with a very hard neutrino 
spectrum. This makes these sources extremely interesting 
for neutrino telescopes. 

The bubbles are parameterized as circular regions of 25° 
radii emitting a uniform flux of neutrinos per steradian. In 
Fig. 2, the relative signal strength expected from each re- 
gion is simulated using the 59-string configuration of the 
IceCube detector. As expected, a part of one bubble that 
dips below the horizon in equatorial coordinates as seen 
by IceCube (upgoing hemisphere) has the largest expected 
signal. Thus the exact sensitivity of IceCube to the Fermi 
bubbles depends on the characterization of the shape of the 
bubbles. However, despite the extensive shape-fittings of 
the gamma emission performed in [18], due to the intrin- 
sic haziness of the structure, the dependence on the fit- 
ting of other emission mechanisms contributing to the ob- 
served gamma-ray structure, and the unknown energy de- 
pendence of the shape at the high-energy region of our in- 
terest, we conclude that using a simple model of two cir- 
cular neutrino-emission regions is adequate for an initial 
LH analysis. The LH analysis will combine data from 
both the 40 and 59-string configurations. Using this cir- 
cular parameterization and assuming a continuous flux of 
3.5 x E 2 dN/dE [1CT 10 TeV cm" 2 s" 1 ], motivated by an 
optimistic gamma-ray flux calculation [20], 164 events are 
expected in the combined 40 and 59-string configuration of 
IceCube. 
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AMANDA Detector 
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Abstract: During the season 2008/2009, IceCube took data as a combined detector with AMANDA embedded into the 
40-string array. With a smaller spacing between the sensors compared to IceCube, AMANDA improved the effective area 
below a few TeV and acted as a first generation low-energy extension of IceCube. The data obtained in this configuration 
is used to search for neutrino sources within the Galaxy. The TeV 7-ray spectra of some potential galactic cosmic-ray 
accelerators show cut-offs in the energy spectrum at energies of a few TeV. In the case of transparent TeV 7-ray sources, 
high-energy neutrinos will follow similar spectra and an improved effective area below a few TeV improves the sensitivity 
for these of sources. 

Several tests including a scan of the galactic plane in the Northern Hemisphere and a dedicated analysis for the Cygnus 
region are presented. In the absence of a significant signal, upper limits are reported. The results provide the most 
restrictive upper limits for the Cygnus region obtained so far. Depending on the assumed energy cut-off, the upper limits 
obtained with this analysis are only a factor of two to three above the expected neutrino flux if all the TeV 7-rays observed 
in the region were of hadronic origin. This implies that during the coming years, IceCube will be able to either detect 
neutrinos from the Cygnus region, or to constrain the nature of the high-energy 7-ray emission in the region, and thus the 
fraction of interacting cosmic rays produced in one of the most active parts of the Galaxy. 
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1 High-Energy Neutrino Production in the 
Galaxy 

One of the primary goals of the IceCube experiment [1] is 
the detection of astrophysical sources of high-energy neu- 
trinos. A neutrino signal uniquely identifies the sites of 
hadron acceleration and interaction and thereby the sites of 
cosmic -ray production. 

If protons are accelerated to sufficiently high energies in 
(galactic) sources, high-energy neutrinos can be produced 
through proton-proton interactions with the ambient gas. If 
such sources are transparent and the 7-ray emission from 
high-energy electrons is small compared to the total 7-ray 
emission, the high-energy neutrino spectra can be inferred 
from the 7-ray spectra [2] . 

Several objects within the Galaxy such as supernova rem- 
nants and pulsar wind nebulae, binary systems and the col- 
lective winds of massive stars might accelerate protons up 
to PeV energies [3]. Even though 7-rays up to several TeV 



have been observed from several of these objects, many 
of the observed 7-ray spectra have energy cut-offs below 
10 TeV and/or their energy spectra are significantly steeper 
than an E~ 2 spectrum as expected from shock acceleration. 
A search for neutrino emission from within the Galaxy thus 
requires an approach that is optimized to retain a high effi- 
ciency for neutrinos with energies below 10 TeV. The pre- 
dicted neutrino flux from galactic sources is very low and 
single point-like sources might elude a detection in the near 
future. In star forming regions, it is however possible that 
several (weak) sources produce an integrated signal strong 
enough for a discovery. In particular the Cygnus region as 
the most active part of the Galaxy in the Northern Hemi- 
sphere is of primary interest to IceCube in this context. 

2 Methods and Targets 

At energies between a few hundred GeV and a few TeV, 
the field of view of the 2008/2009 configuration of Ice- 
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Cube is restricted to the Northern Hemisphere where the 
atmospheric muon background is suppressed by several or- 
ders of magnitude by the shielding provided by the Earth. 
Within this field of view, a search for point-like, steady 
high-energy neutrino sources has been performed. We 
search for a significant excess of neutrinos over the uni- 
form background of atmospheric neutrinos by a maximum 
likelihood ratio hypothesis test, described in [4]. The 
search is performed on a 0.25° grid covering the galac- 
tic plane within the field of view (37.5° < I < 212.5° , 
—5° < b < 5°). Since the angular grid size is smaller than 
the angular resolution of the detector, this search may be 
regarded as an unbinned analysis. The energy term in [4] is 
omitted as it is not relevant for soft spectra sources. 

In addition, the same likelihood ratio hypothesis test is ap- 
plied to six prominent 7-ray sources: the Crab Nebula, LSI 
+61 303, CasA, W51, SS433 and IC443. The sources were 
chosen due to their brightness in 7-rays and/or the presence 
of target material for proton-proton interactions in or near 
the sources. 

To search for high-energy neutrino emission within the 
Cygnus region, a dedicated test based on a 2-point correla- 
tion function has been developed [5]. A search for a spatial 
clustering of events inside a 7°xll° region (72° < I < 83°, 
—3° < b < 4°) around the most active part of the Cygnus 
region is performed. The method is able to take advan- 
tage of extended emission regions and the emission of any 
sources within the region. If applied to the data sample 
used in this work, the discovery flux per point source is 
lower than in a standard search if more than two point 
sources are present within the region. 

The analyzed data set is dominated by atmospheric neu- 
trinos. Any potential astrophysical signal presents only a 
very small contribution in number of events. A data driven 
background estimation can thus be obtained by randomiz- 
ing the arrival directions of the events, compatible with a 
homogeneous spatial distribution. All statistical tests re- 
ported here use this technique for the construction of their 
respective null hypothesis. 

3 The Combined IceCube 40-string and 
AMANDA Detector 

The full IceCube [1] neutrino telescope at the South Pole 
consists of a volume of approximately one cubic kilometer 
of clear Antarctic ice instrumented with light sensors. This 
instrumentation allows to detect muons from charged cur- 
rent interactions of neutrinos. 5160 digital optical modules 
are deployed in the ice along 86 strings that hold 60 optical 
modules each. The detector has been built in several stages 
and new strings have been added each Antarctic summer 
since 2004/2005. 

AMANDA [8] is located at the same site as IceCube and 
consists of 677 optical modules deployed on 19 strings, 
most of them at depths between 1500 and 2000 m. Both 
the vertical and horizontal spacing of the optical modules 



in AMANDA are smaller than in IceCube. This provides 
a lower energy threshold and a higher collection efficiency 
for muons below a few TeV. As AMANDA is the precursor 
experiment to IceCube, many of the techniques employed 
in IceCube were developed and tested in AMANDA [9]. 
AMANDA took data as a stand-alone neutrino telescope 
until December 2006. Since 2007, AMANDA is fully sur- 
rounded by IceCube strings and was integrated into the Ice- 
Cube data acquisition system as a low-energy extension of 
the IceCube detector until 2009 [10]. 

In the combined data taking mode, AMANDA initiates a 
read-out of IceCube whenever a multiplicity trigger con- 
dition in AMANDA is fulfilled. The data collected from 
both parts of the combined detector is then merged into a 
single event and reconstructions can be applied to either the 
full event information or to the IceCube or AMANDA in- 
formation separately. The analysis presented here uses the 
information from both IceCube and AMANDA. 

4 Neutrino Sample 

The targets of this analysis are soft-spectrum sources or 
sources with high-energy cut-offs below the PeV range 
within the Galaxy. The event selection is thus aimed to 
improve the effective area for neutrino energies below this 
scale. This is achieved both through the use of AMANDA 
as an embedded array inside the 40 IceCube strings and 
through an event selection optimized for a larger accep- 
tance of events below 10 TeV compared to the analysis pre- 
sented in [11]. 

The current analysis uses data collected from April 5, 2008 
to May 20, 2009. Both parts of the combined IceCube- 
AMANDA detector operated very stably during this time. 
For IceCube <~ 375 days of data were collected. AMANDA 
was decommissioned before the end of the IceCube 40- 
string data taking period and ^306 days of combined 
IceCube- AMANDA data were collected . The main causes 
for detector downtime were scheduled operations in the 
course of the integration of new strings into the detector. 

The dominant class of recorded events are atmospheric 
muons incident from the atmosphere above the detector. 
The majority of this background is suppressed by a cut on 
the reconstructed direction such that only events from the 
Northern Hemisphere are accepted. Even after this cut, the 
atmospheric muons dominate over the atmospheric neutri- 
nos by several orders of magnitude as a fraction of atmo- 
spheric muons are not well-reconstructed and as such end 
up as up-going muons. In particular coincidences between 
two muons from different air showers can mimic up-going 
events. An event selection is then applied to reduce this 
background by rejecting events with poor reconstruction 
quality and/or events with a high probability to be com- 
posed of two separate particles. An overview of reconstruc- 
tion quality estimators and other event parameters that al- 
low to distinguish signal from background is given in [1 1]. 
A subset of these parameters is used in this work. 
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Figure 1: Energy distribution of atmospheric neutrino 
events in the combined IceCube 40-strings and AMANDA 
point source sample. Events with AMANDA trigger 
("combined events") peak at lower energies than events 
with IceCube trigger. 90% of the events are contained in 
a central interval from 130 GeV to 7.5 TeV. 

The data collected in 2008/2009 contains events that trig- 
gered AMANDA, events that triggered IceCube and events 
that triggered both detectors. Most analyses performed on 
the data sample, such as [11], use only the events that trig- 
gered IceCube. The approach presented here extracts a 
neutrino sample from all three kinds of events. Combined 
IceCube-AMANDA events are selected by different event 
selection criteria than events that only have an IceCube trig- 
ger, as outlined below. 

IceCube events without AMANDA trigger are selected by 
a series of one-dimensional cuts on event quality param- 
eters followed by a multivariate classification based on 
the Neyman-Pearson rule (see for example [12]). The 
probability density functions for five quality parameters 
are generated from atmospheric muon-dominated data as 
background and from atmospheric neutrino simulation and 
combined in the cut. The main cut variable is the likeli- 
hood ratio between the atmospheric neutrino and the atmo- 
spheric muon hypothesis. The distribution of this main cut 
variable is shown in Figure 2 for data and for atmospheric 
neutrino simulation. For combined IceCube-AMANDA 
events, the Neyman-Pearson rule is not applied because a 
series of a series of one-dimensional cuts resulted in a sim- 
ilar performance for these events. 

The energy distribution at the final event selection level 
is shown in Figure 1 . The combined IceCube-AMANDA 
events peak at lower energies. The angular resolution of the 
sample depends on the energy of the events. An unbroken 
E~ 3 power law spectrum has been used to benchmark the 
performance of the analysis. For this very soft spectrum, a 
median angular resolution of 1.2° is achieved. From simu- 
lation of single and double coincident cosmic -ray air show- 
ers with CORSIKA [13], the atmospheric muon contami- 
nation of the cleaned data sample used in this analysis is 
estimated between 2% and 3%. 

19797 neutrino candidates are selected from IceCube and 
AMANDA triggered events. This sample is predominantly 
background atmospheric neutrinos, which cannot be dis- 
tinguished from the cosmic neutrino signal on an event-by- 
event basis 1 . These events are analyzed with the hypoth- 



Figure 2: Distribution of the main cut variable for IceCube- 
triggered events before final cuts are applied. Shown 
are the data in black and atmospheric neutrino simulation 
based on the atmospheric neutrino flux model of [6] in gray. 
The final cut is placed at 1 .0 as result of an optimization of 
the discovery potential [7]. 

esis tests described in the previous section to search for 
spatial clustering of events over the uniform atmospheric 
neutrino background. 

5 Results 

The results of the galactic plane scan are presented in Fig- 
ure 3, from which it is seen that all observations are com- 
patible with the background expectation. The largest clus- 
tering of events was observed at (85.5°, —2, 0°) with a 
(pre-trial) probability to observe an equal or stronger ex- 
cess at this position of 0.0935% due to background fluc- 
tuations only. Accounting for the trials introduced by the 
repetition of the test along the galactic plane, an equiva- 
lent or more significant observation is made in 88.02% of 
randomized data samples. Thus the observed excess in the 
scan is consistent with fluctuations of background. Also 
the results for the six 7-ray sources are compatible with 
the background expectation and preliminary 90% flux up- 
per limits are summarized in Table 1 assuming a power-law 
with a spectral index of 3. The preliminary limits do not 
include the systematic uncertainty of the signal efficiency. 
The strongest preliminary flux limit can be set for Cas A at 
a flux of 5.9 • 10 -11 TeV -1 cm -2 s _1 . The upper limits 
are calculated using the approach of Feldman and Cousins 
[15]. 

With 55 events observed within the box defined around 
the most active part of the Cygnus region compared to 
a background expectation of 60 events, strong flux upper 
limits could be extracted for this region. Assuming an 
E~ 2 - 6 spectrum as was fit to the MILAGRO 7-ray observa- 
tions [16], a preliminary 90% flux upper limit of 3 • 10 -11 
TeV^ 1 cm _2 s _1 (without systematic uncertainties) is ob- 
tained provided the astrophysical signal from the region has 
an exponential energy cutoff only at or above 1000 TeV. 

1 . An exception could be the use of a veto against atmospheric 
neutrinos as proposed in [14] 
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Figure 3: Result of the galactic plane scan using data collected by IceCube 40-strings as a combined detector with 
AMANDA. The significance of the observation at each grid point is expressed by the (pre-trial) pvalue which is shown 
together with the distribution of the events shown as black dots. The most significant excess of events is located at 
(85.5°, -2,0°) with a (pre-trial) pvalue of 0.000935 (-logl0(pvalue)=3.03). The probability to observe a similar or 
stronger excess of events at any point of the galactic Plane is 88.02%. No neutrino sources have been detected. 
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14.15° 
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Table 1: Position, best fit number of source events (ns) 
and preliminary 90% upper limits on the flux of muon 
neutrinos for each of the tested objects in units of 10~ n 
TeV -1 cm~ 2 s _1 . The upper limits are calculated without 
systematic uncertainties under the assumption of an E -3 - 
spectrum. 



6 Outlook 

The installation of the IceCube detector has been com- 
pleted in 2010 and IceCube is now taking data in its fi- 
nal configuration of 86 strings. The collaboration contin- 
ues to search for neutrino sources within and outside of the 
Galaxy. With the substantially larger detector, the sensitiv- 
ity to galactic neutrino sources will improve significantly 
with respect to the analysis presented here. 

In particular IceCube's observations of the Cygnus region 
will enter an interesting regime in the next few years. We 
have shown that the IceCube 40-string/AMANDA limits 
for the Cygnus region are only a factor of two to three 
above the expected flux if all of the 7-rays in the region 
were of hadronic origin. Applying the same test to the data 
obtained with larger configurations of IceCube, it will thus 
be possible to either detect neutrino emission in this region 
or to constrain the hadronic component in the 7-ray emis- 
sion. 

AMANDA has been decommissioned in 2009 and replaced 
by the DeepCore extension of IceCube. The positioning 
of a more densely instrumented volume in the deepest and 
clearest ice around the central IceCube string offers sev- 
eral advantages with respect to AMANDA. In particular, it 



offers the possibility to use veto techniques that allow the 
suppression of downgoing atmospheric muons and might 
open the Southern Hemisphere to neutrino astronomy be- 
low several tens of TeV from the South Pole. 
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Abstract: A time-dependent search for neutrino flares from pre-defined directions in the whole sky is presented. The 
analysis uses a time-clustering algorithm combined with an unbinned maximum likelihood method. This algorithm, 
by including a likelihood signal term describing the contribution of many small clusters of signal-like events, provides 
an effective way for looking for weak neutrino flares over different time-scales. The event selection is optimized to 
maximize the discovery potential for the IceCube 59-string (IC59) detector configuration. Sources are being selected 
based on data in the 0.1 to 100 GeV energy range as provided by the Fermi satellite. Subsequently, periods of interest 
based on electromagnetic data are scanned, over larger time-windows as compared to the duration of the corresponding 
electromagnetic flares. 
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1 Introduction 

Finding neutrino point sources in the sky requires locat- 
ing an excess of events from a particular direction over 
the background of atmospheric neutrinos and muons. Sig- 
nal events might present additional features that distin- 
guish them from background, for example a different en- 
ergy spectrum or time structure. For sources which man- 
ifest large time variations in the emitted electromagnetic 
radiation, the signal-to-noise ratio can be increased by test- 
ing smaller time windows around the flare (time-dependent 
search). Following this idea there are in principle two 
approaches to neutrino time-dependent searches: Trig- 
gered and Untriggered. In the first case we are look- 
ing directly for photon-neutrino correlations using specific 
source lightcurves from Multi-WaveLength (MWL) obser- 
vations [1]. In the second case, followed in this work, 
we perform a generalized search for neutrino flares from 
a pre-selected source list, motivated by (but not directly in 
time coincidence with) MWL observations. This approach 
allows to account for possible time lags between photon 
flares and the associated neutrino flares [2]. 

An untriggered unbinned flare search was first devel- 
oped and applied to IceCube data, using a compact list 
of pre-defined source directions [3]. IceCube is km 3 
scale neutrino detector at the South Pole sensitive to TeV- 
neutrinos [4]. A time-clustering algorithm [3, 5], and an 
unbinned maximum likelihood method [6] are the basis of 



this analysis. Such a method finds the most significant flare 
in a long period. The number of trials coming from all 
combinations of event times is increased, reducing the sig- 
nificance. However, for flares sufficiently shorter than the 
total observation period, the time clustering algorithm is 
more sensitive than a time-integrated analysis. 

In this paper, we propose an extension of the method de- 
scribed in [3]. The proposed algorithm can extract not 
only the most significant flare, but also less significant clus- 
ters of events distributed over several weak flares. These 
weaker flares could have any separation in time and there- 
fore may be very difficult to detect or even undetectable 
with other existing point-sources methods (like [6]). 

2 Multiple flare search algorithm 

A more detailed description of the proposed method and its 
application for multi-flare Monte Carlo simulation can be 
found in [7]. Here we only briefly describe the main steps 
of the proposed algorithm. In order to identify a series of 
weak flares, we first extract all consecutive doublets that 
can be formed out of all signal-like events (Si/Bi > 1) 
over the data taking period ATdata 1 , see Figure 1. 

1. A signal-like event is defined as having Si/Bi > 1, where 
Si and Bi are the signal and background Probability Density 
Function (PDF), respectively, as defined for the time-integrated 
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Figure 1: The basic idea of the time-clustering procedure. 

This step serves to isolate all possible (and smallest) time 
windows (Atj) that compose the signal contribution in the 
tested data sample (the total number being M). We call 
these time windows "data segments". Note, that by using 
only doublets as data segments we do not need any assump- 
tions about the distribution of signal events inside a data 
segment i.e. by definition the time probability p tlmc in the 
signal PDF is uniform in time and is given by p tlmc — — l — . 

Then for each data segment the best estimates of the num- 
ber of signal events h s and source spectral index are 
found by maximizing the one-source likelihood as defined 
in [6]. Then for each data segment the individual value 
of the test statistic TS^Ia*.,- is calculated from the likeli- 
hood ratio of the background-only (null) hypothesis over 
the signal -plus-background hypothesis [6]. All data seg- 
ments are then sorted according to TS^Ia^- In the case 
that real signal events are present, some of these data seg- 
ments will contain the signal events while the rest of them 
are due to background fluctuations. Our aim is to extract 
the optimal (best suited) number of data segments (M opt ) 
which compose the total signal contribution in the overall 
period AT data . 

For this purpose, we used a modification of the single- 
source likelihood function ([6]) by replacing the one- 
source signal term Si by the sum of signal sub-terms over 
m data-segments: 



Si 



Xi Xa 



,Ei,j,Atj) 



E™ i Wi 
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where is a weight which describes the strength (signif- 
icance) of the doublet contained in each data segment. As 
was shown in [7] the test statistic is quite well correlated 
with the true number of injected signal events. Thus we 
take U " I S , \, . 

In order to estimate the optimal number of data segments 
A^opt f° r a given configuration of m segments (starting 
from m = 1) we maximize the modified log(£(n s , -f s , m)) 
with n s and 7 S as free parameters. For a given number m 
the minimization returns the best estimates for the number 
of signal events h s and for the spectral index of the source 
%, and the "global" test statistic is calculated from: 



TS(m) = -2 log 



C(x Sl n s = 0) 
£(x s ,h s ,i s ,m) 



(2) 
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Figure 2: Number of events for a 5<r discovery as a function 
of the total flare duration obtained for one flare search and 
two flares search. See text for more details. 

Then, the optimal number of data segments to be stacked 
(M opt ) is chosen according to the maximum of TS(to). 
The overall significance of the optimal configuration M opt 
can be determined using MC simulations by applying the 
same procedure to a large number of scrambled data sets. 

3 Data Sample 

IceCube 59-string data from May 20, 2009 to May 3 1 , 20 1 
is used. It spans 375 days with an overall effective detector 
uptime of 93% (i.e. 348 days). The whole sky (declina- 
tion range from -85° to 85°) is scanned. Different selec- 
tion criteria, due to different backgrounds, are applied for 
the northern and southern skies, see [8]. After selection 
the data set contains 107 569 events (43 339 events in the 
northern sky and 64 230 events in the southern sky) with a 
median angular resolution of 0.7° . 

4 Expectations for IceCube 59-strings 

Figure 2 shows the performance of the algorithm for a 
source with a E~ 2 energy spectrum located at declina- 
tion (5 = 22°. For a single flare search (solid line with 
open boxes) the number of events for a 5er discovery de- 
creases when we consider flares with shorter duration. As 
an example, for a flare with duration of 28, 10 and 0.1 
days in average about 15, 12 and 7 events, respectively, 
are needed for discovery. Note, that for flares with rela- 
tively short durations (below about 15 days) the number of 
events is smaller compared to a time-integrated analysis, 
see dashed line in Figure 2 labeled: Time Integrated Anal- 
ysis (Space+Energy). In Figure 2 the performance of the 
algorithm for two flare searches is also shown (horizontal 
dashed lines). 



In this case two individual flares with duration At 

,(2) 



(1) 

flare 



and 



Ai .f lare , respectively, are separated in time by a time in- 



method [6]. To calculate the ratio, Si/Bi > 1, only the spatial 
and energy terms in the PDF's are included. 
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terval x. 1 . For two flares separated by x we can see that 
the number of events needed for discovery only slightly de- 
pends on the total flare duration AT{M opt ) and equals the 
case of single flare with duration (AT(M opt ) — Atff are + 

(2) 

^ flare) i- e - x = 0. This is a consequence of the fact 
that the proposed algorithm looks for the total signal in 
the data sample but disregards how these signal events are 
distributed in time. In other words, signal events form in 
time one significant cluster of events (one flare) for a given 
source location, or theses events can be distributed among 
a few (sometimes less significant) flares separated in time. 
Figure 2 shows also that multiple flares have better discov- 
ery potential than that of one flare if the same method is 
used. 

For comparison purposes, in Figure 2 the performance of 
an untriggered time-dependent analysis from [8] is shown. 
In this case calculations are performed using the standard 
unbinned likelihood method with the assumptions, that the 
shape of the flare follows the Gaussian distribution i.e. so 
called Gaussian burst [6] 3 . Comparing our results for 
a single flare with [8] we need about 50% more events 
for discovery in case all events are injected in one single 
flare. This is because our algorithm stacks also background 
fluctuations, and thus leads to a higher 5cr threshold than 
the threshold obtained by a single-source likelihood based 
method. However, if we consider two flares separated in 
time the differences in the number of events strongly de- 
creases and with enough separation in time the multi-flare 
analysis requires fewer events for discovery than standard 
untriggered searches. As an example for two flares 1 day 
long each, if they are separated by more than 20 days, 
the multi-flare search performs better than [8]. A simi- 
lar behavior is observed for individual flares with duration 
^flare ~ 0.005 day, i.e. for time scales of the order of 
minutes. 

In Figure 3 the fluence sensitivity for IceCube 59-string 
data is presented for six representative source directions. 
The fluence depends on the total data period ATdata con- 
sidered, being better for smaller data periods. The fluence 
increases when we consider flares with longer duration. 
The effect is well visible for ATdata = 40 days. 

5 Source selection and results 

The proposed algorithm was applied to selected sources 
which manifest large time variations in the electromagnetic 
flux. Using our multiple search algorithm we do not need 
a precise estimation of the starting time and ending time 
of each flare. As was shown in [7] the algorithm finds all 
signal events in the data period even if signal events are ar- 
ranged as a few clusters separated in time. Thus we only 
need a first guess of the flare central time T m and we set a 
larger time window: T m ± 40 days. This allows to search 
for neutrino flares near a 7, optical, x-ray or infra-red flare 
testing the correlation or anticorrelation in the neutrino- 
gamma emission. 
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Figure 3: Fluence sensitivity from an E~ 2 spectrum neu- 
trino signal plotted versus declination for different obser- 
vation times using IceCube 59-string data. 

In the context of hadronic models predicting high energy 
neutrino emission from objects such as Active Galactic 
Nuclei (AGN), there are several possible scenarios. For 
example in [2] Flat Spectrum Radio Quasars (FSRQ) are 
more promising, as neutrino sources, than BL-Lac objects, 
whereas in [9] the opposite is predicted. The proton blazar 
model [10] predicts that the Low synchrotron peaked BL- 
Lacs (LBL) are more likely to produce a significant neu- 
trino emission than the High synchrotron peaked BL-Lacs 
(HBL). In [9] harder sources are selected as promising 
sources to be detected by IceCube once the prediction of 
the neutrino fluxes, within the assumed pp model, is com- 
bined with the IceCube instrumental response. On the other 
hand in [2] the considered pj model leads to the conclusion 
that FSRQs bright in the GeV range are promising neutrino 
sources without any assumption on the spectral index. In 
order to include these different predictions in this analy- 
sis data from the first Fermi LAT catalog [14] was used to 
select AGNs according to the following criteria: 

• BL-Lacs: Average flux [1 - 100 GeV] > 1 x 10~ 9 
ph cm _2 s _1 AND Spectral index < 2.3 

• FSRQs: Average flux [0.1-1 GeV] > 7 x 10~ 8 ph 
cm _2 s _1 

In addition for both cases we include a variability index 
cut (V > 23.21) to select sources that are more likely to 
exhibit flaring periods [14]. For the selected sources, infor- 
mation about flaring states in different wavelengths, during 
the period of the 59-string configuration of IceCube was 
collected. The selected sources and periods are presented 
in Table 1 . 



2. In this case the total flare duration is defined as 



At 



(i) 



f 2") 

+ x + Atj- lare , so for example config- 



AT(M ov t) = ^ flare 

uration 1 + x + 1 corresponds to two one day flare each separated 
by time interval x ranging from up to 27 days. 

3. In this work simulated signal events are injected accord- 
ing to a uniform time distribution, while in [8] a Gaussian dis- 
tribution with different standard deviation a w was considered. 
To make comparison, the corresponding width at half maximum 
FWHM = 2\/2 \u2a w is calculated for a Gaussian flare, and 
then these results are compared with our calculations. 
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Table 1 : Results for pre-defined variable astrophysical source candidates using the multi-flare algorithm. 
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* T m is the midpoint of the flare time interval reported in the alert (Atel) or in the corresponding reference when available. 
AT(M pt) is the flare duration calculated for the optimal configuration of M op t data segments. 
The fluence upper limit is calculated by integrating d$/dE x E over the 90% energy range and AT(M opt ), assuming a neutrino 

energy spectrum of E~ 2 . 



The proposed method was applied to the selected source 
candidates.No significant excess above the atmospheric 
background is found, therefore upper limits on the neutrino 
fluence were calculated. The results are presented in Ta- 
ble 1 . The highest fluctuation observed corresponds to 3C 
454.3 with a p-value of 8% (not including the trial factors 
due to looking at several sources). The limits for IC59 are 
on average about 50% better than for IceCube 40-strings 
data [15]. 

6 Summary 

We presented a method to search for neutrino flares from 
point sources without an a-priori assumed time structure. 
The method considers only data segments which contain 
signal-like doublets, and uses a test-statistic term as their 
weights in a stacking-like calculation for the global max- 
imum likelihood. For flares sufficiently shorter than the 
total observation period, the method is more sensitive than 
a time-integrated analysis and in some cases is also more 
sensitive than single flare searches already in use (like [8]). 
IceCube 59-string data was analyzed using the proposed 
method looking for neutrino multi-flares with no a-priori 
assumption on the time structure of the signal. A list 
of promising source candidates was selected based on 



different hadronic models. Since no deviation from the 
background-only hypothesis was found, upper limits on the 
neutrino fluence from these sources were derived. 
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Abstract: Transient neutrino sources such as Gamma-Ray Bursts (GRBs) and Supernovae (SNe) are hypothesized to 
emit bursts of high-energy neutrinos on a time-scale of < 100 s. To increase the sensitivity to detect those neutrinos 
and identify their sources, an optical follow-up program for neutrinos detected with the IceCube observatory has 
been implemented. If a neutrino multiplet, i.e. two or more neutrinos from the same direction within 100 s, is found 
by IceCube a trigger is sent to the Robotic Optical Transient Search Experiment, ROTSE. The 4 ROTSE telescopes 
immediately observe the corresponding region in the sky in order to detect an optical counterpart to the neutrino events 
of IceCube. Data from the first year of operation of the optical follow-up program have been searched for a signal 
from supernovae. No statistically significant excess in the rate of neutrino multiplets has been observed and further no 
coincidence with an optical counterpart was found during the first year of data taking. This allows us to restrict current 
models predicting a high-energy neutrino flux from soft jets in core-collapse SNe. For the first time a stringent limit on 
the hadronic jet production in core-collapse SNe is derived. 
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1 Introduction 

High-energy astrophysical neutrinos are produced in pro- 
ton interactions of charged cosmic rays with ambient pho- 
ton or baryonic fields (for reviews see [1]). Acceleration 
of protons to very high energies takes place in astrophys- 
ical shocks. Neutrinos escape the acceleration region and 
propagate through space without interaction, while protons 
are deflected in magnetic fields and no longer point back 
to their source. Unlike gamma-rays, neutrinos are solely 
produced in hadronic processes and could therefore re- 
veal the sources of the highest energy charged cosmic rays. 
Gamma-ray bursts could provide the environment and the 
required energy to explain the production of the highest en- 
ergy cosmic -rays and hence are a plausible candidate. Re- 
cent observations imply a common physical origin of long 
GRBs and core-collapse supernovae (CCSNe): a massive 
stellar explosion (see [2] for a review). According to the 
collapsar model [3], long GRBs (duration > 2 s) have their 
origin in the collapse of a massive, rapidly rotating star into 
a black hole surrounded by an accretion disk. Relativistic 
jets with Lorentz boost factors of 100-1000 form along the 
stellar axis. This GRB-SN connection gives rise to the idea 
that GRBs and SNe might have the jet signature in com- 
mon and a certain fraction of core-collapse SNe might host 
soft relativistic jets. SN jets are suggested to be equally en- 



ergetic and more baryon-rich, hence they are only mildly 
relativistic. Such soft relativistic jets would become stalled 
in the outer layers of the progenitor star, leading to essen- 
tially full absorption of the electromagnetic radiation emit- 
ted by the jet and at the same time an efficient production 
of high-energy neutrinos [4, 5]. This motivates a search 
for neutrino emission, as neutrinos would be able to escape 
from within the star. 

The IceCube neutrino detector, located at the geographic 
South Pole, is built to detect high-energy astrophysical neu- 
trinos [6]. So far GRB neutrino searches have been per- 
formed offline on AMANDA [7] and IceCube [8] data, trig- 
gered by gamma-ray satellite detections. Furthermore, a 
dedicated search for a neutrino signal in coincidence with 
the observed X-ray flash of SN 2008D has been conducted 
by IceCube [9] in order to test the soft jet scenario for CC- 
SNe. Neither the GRB nor the SN neutrino searches led to 
a detection yet, but set upper limits to the possible neutrino 
flux. 

Early SN detections, as in the case of SN 2008D, are very 
rare since X-ray telescopes have a limited field of view. 
However, neutrino telescopes cover half of the sky at any 
time. If neutrinos produced in soft relativistic SN jets are 
detected in real time, they can be used to trigger follow-up 
observations [10]. This is realized with the optical follow- 
up program presented here. Complementary to the offline 
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searches, the optical follow-up program is an online search 
independent of satellite detections. It is sensitive to tran- 
sient objects, which are either gamma-dark or missed by 
gamma-ray satellites. In addition to a gain in significance, 
the optical observations may allow to identify the transient 
neutrino source, be it a SN, GRB or any other transient phe- 
nomenon producing an optical signal. Hence it enables us 
to test the plausible hypothesis of a soft relativistic SN jet 
and sheds light on the connection between GRBs, SNe and 
relativistic jets. 

In order to implement the optical follow-up program an on- 
line neutrino event selection was developed at the neutrino 
detector IceCube. The data are processed online by a com- 
puter farm at the South Pole. A multiplicity trigger selects 
neutrino burst candidates and the directional information is 
transferred to the four ROTSE telescopes, which start the 
follow-up immediately and continue observations for sev- 
eral days. The obtained optical data are analyzed in order 
to search for an optical supernova counterpart. 

2 IceCube 

The IceCube neutrino telescope has been under construc- 
tion at the geographic South Pole since 2004 and was com- 
pleted in the Antarctic summer of 2010/11. It is capa- 
ble of detecting high energy neutrinos with energies above 
100 GeV and is most sensitive to muon neutrinos within the 
energy range from TeV to PeV. High-energy muon neutri- 
nos undergoing charged current interactions in the ice or 
the underlying rock produce muons in neutrino-nucleon in- 
teractions. The muon travels in a direction close to that 
of the neutrino and emits Cherenkov light. The deep ultra 
clear Antarctic ice is instrumented with light sensors thus 
forming a Cherenkov particle detector. After its completion 
it comprises a volume of 1km 3 with 5160 digital optical 
modules (DOMs) attached to 86 vertical strings at a depth 
of 1450 m to 2450 m [6]. Each DOM consists of a 25 cm 
diameter Hamamatsu photomultiplier tube (PMT) and sup- 
porting hardware inside a glass pressure sphere. Here we 
present the analysis of the data taken from the start of the 
follow-up program on 2008/12/16 to 2009/12/31. Initially 
40 IceCube strings were taking data. In May 2009 an ad- 
ditional 19 strings were included. This corresponds to an 
uptime of 121 days with 40 and 186.4days with 59 strings. 
In the following the deployment stages will be referred to 
as IC40 and IC59. 

2.1 Online System 

In order to rapidly trigger optical telescopes the first online 
analysis of high-energy neutrinos detected by IceCube was 
developed and implemented. Unlike in the offline analyses, 
which are performed on an entire dataset (usually ~ 1 year 
of data) with time consuming reconstructions on a large 
computer cluster, the data are processed online by a com- 
puter cluster at the South Pole. The processing includes 
event reconstruction and basic event selection. The first 



year of data presented here was taken with a latency of 6- 
8 h. With the start of operations with 79 strings the pro- 
cessing was upgraded reducing the latency to a few min- 
utes. After the parallel processing the data arrive on a dedi- 
cated machine (analysis client), where a sophisticated event 
selection is applied based on the reconstructed event pa- 
rameters. A multiplicity trigger selects neutrino burst can- 
didates (see section 2.2). No further reconstruction algo- 
rithms need to be applied at the analysis client allowing a 
very fast filtering of the events (<C 1 s). The directional 
information is transferred to Madison, Wisconsin, via the 
Iridium satellite network within about 10 s. From there the 
message is forwarded to the four ROTSE telescopes via the 
internet through a TCP-socket connection for immediate 
follow-up observations. The stability and performance of 
the online system is constantly monitored in order to allow 
a fast discovery of problems. To achieve this, test alerts are 
produced at a much higher rate (~ 100 test alerts per day 
compared to 25 real alerts per year) by the same pipeline 
and are also send to the North. Their rate and delay time 
distributions are monitored using an automatically gener- 
ated web page. 

2.2 Neutrino Event Selection 

The background in a search for muon-neutrinos of astro- 
physical origin can be divided into two classes. One con- 
sists of atmospheric muons, created in cosmic ray air show- 
ers, entering the detector from above. The other is given 
by atmospheric neutrinos which have their origin in meson 
decays in cosmic ray air showers. The expected neutrino 
signal according to the soft jet SN model can be calculated 
as a function of two model parameters: the boost Lorentz 
factor T and the jet energy Ej et [9]. Signal events are simu- 
lated following the predicted neutrino flux spectrum in or- 
der to develop and optimize selection criteria to distinguish 
signal and background events. Restricting the search to 
the Northern hemisphere and imposing requirements on the 
event reconstruction quality (e.g. the number of hits with 
small time residual or the likelihood of the reconstruction) 
allows a suppression of the mis -reconstructed muon back- 
ground. To suppress the background of atmospheric neutri- 
nos, which we cannot distinguish from the soft SN neutrino 
spectrum, we require the detection of at least two events 
within 100 s and an angular difference between their two 
reconstructed directions of A\P < 4°. The choice of the 
time window size is motivated by the jet penetration time. 
The observed gamma-ray emission from long GRBs has 
a typical length of 50 s, which roughly corresponds to the 
time for a highly relativistic jet to penetrate the stellar en- 
velope. The angular window A^> is determined by the an- 
gular resolution of IceCube and was optimized along with 
the other selection parameters. The final set of selection 
cuts has been optimized in order to reach a multiplet rate 
of <~ 25 per year corresponding to the maximal number of 
alerts accepted by ROTSE. The final data stream consists of 
37% (70%) atmospheric neutrinos for IC40 (IC59). Com- 
bining the neutrino measurement with the optical measure- 
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ment allows the cuts to be relaxed yielding a larger back- 
ground contamination and at the same time a higher signal 
passing rate. A doublet is not significant by itself, but may 
become significant when the optical information is added. 
Each multiplet is forwarded to the ROTSE telescopes. The 
doublet direction is calculated as a weighted mean from 
the single reconstructed directions comprising the multi- 
plet. The single events are weighted with l/cx 2 , where a 
is the reconstruction error estimated by the paraboloid fit, 
which fits a paraboloid to the likelihood landscape around 
the minimum defined by the best fit. The resolution of the 
doublet direction is ~ 0.8°. 

3 Search for Optical Counterparts 

The IceCube multiplet alerts are forwarded to the robotic 
optical transient search experiment (ROTSE), which con- 
sists of four identical telescopes located in Australia, Texas, 
Namibia and Turkey [11]. The telescopes stand out because 
of their large field of view (FoV) of 1.85° x 1.85° and a 
rapid response with a typical telescope slew time of 4 sec to 
move the telescope from the standby position to the desired 
position. The telescopes have a parabolic primary mirror 
with a diameter of 45 cm. To be sensitive to weak sources 
no bandwidth filter is used. ROTSE is most sensitive in the 
R-band (~ 650 nm). The wide field of view is imaged onto 
a back-illuminated thinned CCD with 2048 x 2048 13.5 /im 
pixels. For a 60 sec exposure at optimal conditions the lim- 
iting magnitude is around mji ss 18.5, which is well suited 
for a study of GRB afterglows during the first hour or more 
and SN light curves with peak magnitude < 16. The cor- 
responding FWHM (full width at half maximum) of the 
stellar images is less than 2.5 pixels (8.1 arcseconds). Ob- 
servations are scheduled in a queue and are processed in 
the order of their assigned priority. IceCube triggers have 
second highest priority after GRB follow-ups triggered by 
the GRB Coordinate Network (GCN). 
Once an IceCube alert is received by one of the tele- 
scopes, the corresponding region of the night sky will be 
observed within seconds. A predefined observation pro- 
gram is started: The prompt observation includes thirty ex- 
posures of 60 seconds length. Follow-up observations are 
performed for 14 nights. This was extended on 2009/10/27 
to 24 nights, with daily observations for 12 nights and then 
observations during every second night up to day 24 after 
the trigger was received. Eight images with 60 seconds ex- 
posure time are taken per night. The prompt observation is 
motivated by the typical rapidly decaying light curve of a 
GRB afterglow, while the follow-up observation of 14 (or 
24) nights permits the identification of a rising SN light 
curve. In the initial phase with IC40 and IC59, the online 
processing latency of several hours made the search for an 
optical GRB afterglow unfeasible. We therefore focus on 
the SN light curve detection in the ROTSE data. 
Image correction and calibration are performed at the tele- 
scope sites. The images of each night are combined in order 
to obtain a deeper image. A reference image is subtracted 



from each combined image using the algorithm developed 
by [12]. As deep images are usually not available for the 
positions we would like to observe, we initially choose the 
deepest image of our observing sequence as the reference 
image. In 40% of the alerts we took another deep image 
roughly one year later. Both SN light curves and GRB af- 
terglows would have faded after a few weeks, and would 
not be present in the newly taken reference image. 
All extracted objects found in the subtracted images are 
candidates for variable sources. However, bad image qual- 
ity, failed image convolution, bad pixels and other effects 
frequently cause artifacts in the subtraction process, requir- 
ing further selection of the candidates. A candidate identi- 
fication algorithm including a boosted decision tree is ap- 
plied to classify sources according to geometrical and vari- 
ability criteria. The final candidates are summarized on a 
web page and are inspected visually by several trained per- 
sons, who have to classifly the candidate as a SN, a variable 
star or a subtraction artifact. SN candidate identification by 
the human eye works well as shown in the galaxy zoo SN 
project [13]. The visual scanning was performed by three 
individual persons to ensure no good candidate was missed 
and to avoid false positives. 

4 Results 

This paper presents the results from the analysis of data 
taking in the period of 2008/12/16 to 2009/12/31. Ta- 
ble 1 shows the number of detected and expected doublets 
and triplets for the IC40 and the IC59 datasets as well as 
the number of detected and expected optical SN counter- 
parts. The IceCube expectation based on a background 
only hypothesis was obtained from scrambled datasets. 
To correctly incorporate detector asymmetries, seasonal 
variations and up-time gaps we used the entire IC40 and 
IC59 datasets and exchanged the event directions ran- 
domly while keeping the event times fixed. The number 
of doublets shows a small excess, which corresponds to 
a 2. 1 a effect and is thus not statistically significant. The 
expected number of randomly coincident SN detections, 
iV s ^ = 0.074, is based on an assumed core-collapse SN 
rate of 1 per year within a sphere with radius lOMpc, i.e. 
2.4 • 10~ 4 y _1 Mpc~ 3 , and a Gaussian absolute magnitude 
distribution with mean of —18 mag and standard deviation 
of 1 mag for CCSN [14]. In total 31 alerts were forwarded 
to the ROTSE telescopes. Five could not be observed be- 
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Figure 1: Limits on the choked jet SN model [5] for dif- 
ferent boost Lorentz factors T as a function of the rate of 
SNe with jets p and the jet energy E- ]ct (colored regions are 
excluded at 90% CL). Horizontal dashed lines indicate a 
fraction of SNe with jets of 100%, 10% or 1%. 

cause they were too close to the sun. For two alerts no good 
data could be collected. Seven alerts were discarded be- 
cause the corresponding fields were too close to the galac- 
tic plane and hence too crowded. Thus 17 good optical 
datasets remained for the analysis. The data were pro- 
cessed as described above. No optical SN counterpart was 
found in the data. 

We obtain the confidence level for different combinations 
of SN model parameters [5] by using a pre-defined test 
statistic based on a likelihood function. The limit is cal- 
culated for the jet boost Lorentz factors T = 6, 8, 10 as 
a function of the rate of SNe with jets p and the jet en- 
ergy Ej et . Systematic errors related to the simulated neu- 
trino sensitivity and the SN sensitivity are included in the 
limit calculation. The 90% confidence regions for each T- 
value are displayed in the £j et -/9-plane in figure 1 (colored 
regions are excluded at 90% CL). Including the optical in- 
formation into the limit calculation improved the limit and 
allows tests of 5-25% smaller CCSN rates. The largest 
improvement is obtained for small jet energies and large 
CCSN rates. The most stringent limit can be set for high 
T-factors. Less than 4.2% of all SNe have a jet with T = 10 
and a typical jet energy of E- ]ct — 3 • 10 51 erg. This is the 
first limit on CCSN jets using neutrino information. 

5 Summary and Outlook 

The optical follow-up program of IceCube neutrino mul- 
tiplets realized by the four ROTSE telescopes proves the 
feasibility of the program. The technical challenge of an- 
alyzing neutrino data in real time at the remote location of 
the South Pole and triggering optical telescopes has been 
solved. First meaningful limits to the SN slow-jet hypoth- 
esis could be derived already after the first year of oper- 
ation. Especially in cases of high boost Lorentz factors 
of r = 10 stringent limits on the soft jet SN model are 



obtained. Soderberg et al. [15] obtain an estimate on the 
fraction of SNe harboring a central engine from a radio sur- 
vey of type Ibc SNe. They conclude that the rate is about 
1%, consistent with the inferred rate of nearby GRBs. Our 
approach is completely independent and for the first time 
directly tests hadronic acceleration in CCSN, while the ra- 
dio counterpart is sensitive to leptonic acceleration. 
The volume of the IceCube detector has now increased 
to a cubic kilometer yielding a larger sensitivity to high- 
energy neutrinos. In addition the acquired uptime is grow- 
ing continuously. The delay of processing neutrino data at 
the South Pole has been reduced significantly from several 
hours to a few minutes. This results in the possibility of a 
very fast follow-up and allows the detection of GRB after- 
glows, which fade rapidly below the telescope's detection 
threshold. 

Because of the successful operation of the optical follow- 
up program with ROTSE, the program was extended in Au- 
gust 2010 to the Palomar Transient Factory (PTF) [16], 
which will provide deeper images and a fast processing 
pipeline including a spectroscopic follow-up of interesting 
SN candidates. Furthermore an X-ray follow-up by the 
Swift satellite of the most significant multiplets has been 
set up and started operations in February 201 1 [17]. 
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Abstract: The search for neutrinos of astrophysical origin is among the primary goals of the IceCube neutrino telescope. 
Point source candidates include galactic objects such as supernova remnants (SNRs) as well as extragalactic objects 
such as Active Galactic Nuclei (AGN), Supernovae (SNe) and Gamma-Ray Bursts (GRBs). To increase the sensitivity 
of the search for high-energy neutrinos from SNe and especially GRBs an X-ray follow-up with the Swift satellite has 
been developed. Triggered by interesting IceCube events the satellite will be repointed aiming for the detection of a 
transient X-ray counterpart, e.g. an X-ray GRB afterglow. In addition to typical GRBs the program is sensitive to SN 
shock breakouts, slightly off-axis GRBs and orphan GRB afterglows. The online event selection in IceCube as well as 
the X-ray observation strategy will be presented. 
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1 Introduction 

1.1 Overview of the program 

The X-ray Follow-Up with Swift is a multimessanger ap- 
proach developed by the Swift and the IceCube collabora- 
tion to detect GRBs. It uses the IceCube neutrino telescope 
at the South Pole to look for signal like neutrino-multiplets 
(i.e. at least two neutrinos from the same direction within 
100 seconds) to trigger for follow-up observations with the 
Swift satellite in the X-ray band. The implementation of 
the program makes use of the existing neutrino event selec- 
tion of the Optical Follow-Up Program [1] (OFUP) at the 
South Pole. Neutrino multiplets are found online in quasi- 
real time with a typical latency of about 5 minutes. This 
low latency opens the possibility to search for fast decaying 
X-ray afterglows from GammaRaz Bursts (GRBs). Addi- 
tional latency is expected on the Swift side, due to commu- 
nication constraints with the spacecraft, the orbital position 
of Swift, and human-in-the-loop requirements for space- 
craft commanding. Depending on the visibility of Swift to 
a ground relay station, the additional delay will be between 
30 minutes to 4 hours. It is worth noting that the typical X- 
ray afterglow associated with long GRBs is visible to Swift 
for days and sometimes weeks. 



1.2 Scientific Motivation 

Astrophysical neutrino bursts are not necessarily always 
accompanied by an observable prompt electromagnetic 
flux but can be detected using neutrinos. Such objects 
could be GRBs whose narrow jets don't point directly to- 
wards earth [2] or choked GRBs [3] for which the jet may 
fail to penetrate the stellar envelope. Despite the lack of 
prompt 7-rays, the source is likely to be visible in X-rays 
from shock breakout after 10 3 s, and exhibit an optical 
counterpart similar to that seen in core collapse supernovae 
[4]. While no firm estimate exists, the number of these 
dark bursts might be ten times larger than 7-bright bursts. 
Therefore, the search for transient neutrino sources can to 
play an important part in the search for (7-dark) GRBs, and 
may provide insights e.g. into the origin of the high en- 
ergy cosmic rays. While an optical follow-up can be con- 
ducted by ground based telescopes (e.g. ROTSE [5]), Swift 
is uniquely capable of rapid follow up with X-ray observa- 
tions. 

NASA's Swift Explorer Mission is an ideal tool for study- 
ing the electromagnetic radiation from violent astrophysi- 
cal events, such as GRBs. Three telescopes are supported 
by the Swift platform. A wide field of view instrument, 
the Burst Alert Telescope (BAT), monitors for the prompt 
7-rays from a GRB. In response to a burst, Swift will slew 
into position to image the BAT error region with the X-Ray 
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Telescope (XRT) and Ultra Violet and Optical Telescope 
(UVOT). 

2 Alert chain 

Swift is in high demand amongst the scientific community. 
Hence only limited observing time is available for the Ice- 
Cube follow-up program. Using the setup of the OFUP, 
one obtains 25 multiple ts triggers per year with IceCube. 
Most of these are due to atmospheric neutrinos and some 
atmospheric muons. The number of extensive follow-up 
observations can be reduced in two steps to one per year. 
The number of IceCube alerts can be decreased to approxi- 
mately 7/year with very little loss in signal efficiency. This 
is achieved by making a likelihood method, as described in 
section 2.1. The second step involves a quick evaluation of 
the first available X-ray data (section 3.1), before triggering 
the extensive follow-up program. 

2.1 IceCube Trigger Selection 

For the optical follow-up program, the singlet data rate is 
reduced by the optical follow-up online level 3 filter [5] to 
R s rts 2mHz, almost reaching a pure (atmospheric) neu- 
trino sample. Using this data sample, multiplets are se- 
lected if they arrive within At = 100s and from the same 
direction within the reconstruction uncertainty of A* = 
3.5°. The trigger conditions reduce the detected number 
of coincident neutrinos from the isotropic background of 
atmospheric neutrinos to 25 false positives per year. 

To maximize the discovery potential of the IceCube-Swift 
program a new test statistic was developed to test for the 
possibility that a neutrino doublet is of astrophysical ori- 
gin and lies within Swift's field of view. As the derivation 
is beyond the scope of this paper, it is only described and 
motivated here. We begin with the following definitions: 

<T q = <7 1 + cr 2 
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where ui/2 are the reconstruction uncertainties of the par- 
ticipating neutrinos that arrive from the (reconstructed) di- 
rections V\ j2 with an angular difference of if). Assuming a 
circular follow-up region, the test statistic 
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tends to small values for signal-like doublets and larger val- 
ues for background-like events. It takes various effects into 
account. 

• The first two terms act together. While the first term 
favors events with a small angle ip, indicating neu- 
trinos from the same direction and possibly source, 



it also introduces a punishment for small reconstruc- 
tion uncertainties. The qualitative explanation is that 
two neutrinos for which the error regions do not over- 
lap are more likely background than signal. As a 
consequence the first term tends also to small values 
for large combined reconstruction uncertainties a q . 
The second term counteracts this effect, introducing 
a punishment for large uncertainties. Thus, the two 
first terms favor well reconstructed events from the 
same direction. 

• The third term introduces the tiled Swift field of view 
with a radius of Q A w 0.5°. It favors those events 
with small errors for which the reconstructed dou- 
blet direction matches well with the true direction, 
thus minimizing the possibility of observing a region 
of space during a follow-up which does not include 
the actual source within the FoV and supporting the 
first two terms in selecting well reconstructed events. 
Obviously, the Swift FoV need to be evaluated before 
a cut decision can be done. 

• The time difference At between two neutrinos is 
considered in the fourth term. Normalized to the 
100 s time window of the trigger, small values are 
reached for small time differences assuming they are 
an indicator for a neutrino bundle of an astrophysical 
source. 

2.2 Swift Follow-Up 

IceCube provides a median position resolution for selected 
events of less than one degree. However, the XRT field 
of view is only 0.4° in diameter, which will cover only a 
fraction of roughly 20% of the IceCube space angle distri- 
bution. Due to the limited coverage by the XRT field of 
view, we are forced to tile the follow-up region with seven 
pointings of Swift, thus creating an artificially larger field 
of view. In this analysis we treat the regions where two 
fields overlap in an identical way to the other regions (i.e. 
we do not count them twice). The discovery potential for 
bright sources is not expected to improve for longer expo- 
sures and we choose to optimize the follow-up program for 
bright sources that are easily discovered in less than 1 ks. 
With this, a wider field is preferred over a deeper exposure. 

Figure 1 shows the acceptance of an astrophysical source 
with an E~ 2 spectrum by the test statistic relative to the 
optical follow-up filter as a function of the estimated num- 
ber of doublet triggers per year. Starting with loose cuts on 
the test statistic and 25 false positives, one can see the drop 
in the number of false positives as well as signal efficiency 
for tighter cuts reducing the background to zero doublets 
per year. The maximum acceptance equals 63% which is 
the fraction of IceCube's point spread function covered by 
the tiled follow-up region. 

We decided to aim for seven observable alerts per year. 
This reduces the background by 72% while, in the case of 
an E~ 2 spectrum, there is only a loss of 10% (14%) for 
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Figure 1: The acceptance of an astrophysical source with 
an E~ 2 spectrum of the test statistic relative to the opti- 
cal follow-up filter as a function of the estimated number 
of doublet triggers per year. The effect of the Swift field 
of view is included, only accepting events within the FoV 
and reducing the acceptance to a maximum of 63%. It is 
displayed for four different cases assuming a signal with 
a time difference between the arriving signal neutrinos of 
At. The aim is to reduce the background to seven alerts per 
year. 



events arriving with a time difference of 50 (100) seconds 
(figure 1) relative to the maximum acceptance. Almost all 
doublets with a short time difference will pass the cut. 

3 Expected Results 

Swift orbits the Earth every 96 minutes, with the Ice- 
Cube trigger region becoming visible each time for approx- 
imately 2 ks. The spacecraft will be commanded to auto- 
matically observe the seven fields as soon as they rise above 
the Earth limb, providing approximately 285 s of observa- 
tion time each per orbit. The observations will be repeated 
every orbit until a total of approximately 2 ks per tile is 
achieved, typically taking between 12 hours and a full day 
to complete. 

In this way, any X-ray sources visible to Swift will be 
observed multiple times, generating a light curve. These 
products provide critical information for interpreting the 
nature of the source, possibly identifying a GRB via a typ- 
ical GRB afterglow. 

This might prove vital in the process of distinguishing be- 
tween background sources accidentally found within the 
FoV and interesting transient objects. Two different sig- 
nificance tests are proposed. The level 1 test (section 3.1) 
provides a relatively quick test on the first day of data col- 
lected by Swift, to decide whether to initiate a multi-day 
follow-up program. The level 2 test (section 3.2), made on 
the full data set, provides a threshold for claiming a joint 
Swift-IceCube discovery of an X-ray afterglow in coinci- 



dence with an astrophysical neutrino source. The level 1 
test is presented in its final form here, while the level 2 test 
is still under development. 

3.1 Level-1 Significance Test 

After several orbits, Swift may have detected one or more 
X-ray sources, with a position uncertainty that is typically 
on the order of a few arcsec, limited by the Swift XRT 
point spread function. The source position and an initial 
measurement of the flux will allow for a preliminary test 
of whether it stands out from the expected background of 
X-ray sources. The expected number of background X-ray 
sources depends strongly on the flux threshold. We will 
consider an X-ray source to have passed the level 1 signif- 
icance test if it satisfies any of the following criteria for an 
(extra) galactic search. 

A Uncatalogued Sources: The level 1 source is not 
in proximity to a catalogued X-ray object (i.e. not 
within 3 a of the combined Swift and catalog po- 
sition uncertainty), is brighter than a flux threshold 
Sa = (5-10~ 12 )l-10~ 10 erg/(cm 2 s), and occurs in 
a region of the sky where the ROSAT Bright Source 
Catalog would have observed it had it been in its cur- 
rent state when surveyed. 

B Variable Sources: The source is brighter than a flux 
threshold S B = (5 ■ 10~ 13 )1 ■ lO" 11 erg/(cm 2 s) 
and exhibits significant variability across the first day 
of Swift data, with the p-value of a fit to a flat light 
curve being lower than some critical value P. 

C Active Catalogued Sources: The level 1 source is 
within 3tr uncertainty of the position of a catalogued 
object but the new measurement is Mx brighter than 
it appears in the catalog. 

D Poorly Catalogued Sources: A source lies outside 
of the region covered by the ROSAT Bright Source 
Catalog (due to the low exposure time of ROSAT in 
that region) but is observed to be M x brighter than 
the threshold set by ROSAT for that region. 

It is noted that there are significantly more serendipitous 
X-ray backgrounds in proximity to the galactic plane. As 
such, the analysis is carried out with different thresholds 
depending on galactic latitude (GL). Specifically, a higher 
threshold is used if \GL\ < 20°. 

Our intended goal is that the level 1 test will allow no more 
than one false positive per year. This is achieved by ad- 
justing the above parameters (Sa,Sb,P,N,M) to keep 
the number of serendipitous X-ray sources to fewer than 
Ns < 0.18 deg~ 2 .If an X-ray source passes the level 1 
test, then additional Swift data is accumulated over the fol- 
lowing « 1 week, to determine a light curve and spectrum 
for the source. The full data set is then utilized for the level 
2 test. 
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The level 1 tests will first be applied once two observations 
have been made on each tile (i.e. after two Swift orbits). 
The analysis will continue on the accumulating data until 
the level 2 observations are triggered or a total of 2 ks have 
been observed for each tile, whichever is sooner. If a total 
of 2 ks has been accumulated for each tile without trigger- 
ing the level 2 observations, then the fields will be declared 
uninteresting and no further observations or analysis will 
be carried out. 

3.2 Level-2 Significance Test 

Should an X-ray source be discovered that passes the level 
1 test described above, a dedicated observing program will 
be initiated for that source. The tiled observations will be 
discontinued and Swift will take up a pointed observing 
mode, with the source at the center of a single XRT field 
of view (0.4° diameter). The level 2 test will determine 
the significance of all data accumulated over w 1 week, in 
conjunction with the IceCube trigger data. Backgrounds 
will be significantly reduced from level 1 to level 2 by ex- 
amining the larger data set, looking at additional features 
for transient behavior like the slope of the light curve, the 
shape of the X-ray spectrum and the source fading into ob- 
scurity. 

It is anticipated that, for a given class of transient X-ray 
sources (GRB afterglow, AGN activity, etc), we will be 
limited by an irreducible background of similar events. A 
first estimate, based on the rate of BAT-triggered GRBs and 
the average light curve behavior, predicts that serendipi- 
tous GRB afterglows will be discovered only once per 3000 
years with the Swift-IceCube program. 

However, significant effort will be required to put a final 
limit on these chance events, requiring both an analysis of 
previously observed fields and a careful study of theoreti- 
cal constraints to place limits on the number of untriggered 
GRBs like failed or j-dark GRBs, off-axis or orphaned 
GRBs or ordinary GRBs that were not in the field of view 
of a telescope during the prompt outburst, but nonetheless 
produce a serendipitous X-ray afterglow. 

Numerous studies have placed limits on the number of un- 
triggered GRBs, typically of order 100 times the rate of 
regular GRBs [6, 7, 8]. This would place a limit on the 
Swift-IceCube level-2 false positive rate of once per 30 
years. 



Considering the limitations in our Swift search to brighter 
GRBs, we use a conservative probability of 40% instead 
of the 90% detection efficiency of Swift for all X-ray af- 
terglows. The OFUP has an detection efficiency of 40% 
for single neutrino events. Combining these two numbers 
with the fact that about 20% of the sky are not observable 
by Swift due to the sun, 13% of all GRB neutrino events 
could be identified as such [9]. Assuming the prediction of 
10 GRB neutrinos in IceCube per year and hemisphere and 
requiring a doublet detection (introducing a reduction by a 
factor of 20 [9]) one obtains a detection rate of 0.064 GRBs 
per year of which 0.04 GRBs would lie within the FoV of 
Swift (PfoV = 63%). However, according to [6] there 
could be up to 5 times more 7-dark GRBs with an X-ray 
flux above our level 1 threshold of 5 • 10~ 12 erg/(cm 2 s), 
increasing the number of expected detections possibly to 
about 0.2 GRBs per year. 

A future extension triggering on high energy single neu- 
trino events could overcome the reduction factor of 20 for 
doublet observations and add a signal which would be an 
order of magnitude higher, (these numbers might be up- 
dated and change slightly in the next week) 

References 

[1] A. Franckowiak for the IceCube collaboration, ICRC 
proceedings, 201 1 

[2] R. Yamazaki, K. Ioka, T. Nakamura, K. Toma , Adv. 
Sp. Res., 2006, 38:1299 

[3] P. Meszaros, Rept.Prog.Phys., 2006, 69: 2259-2322 

[4] S. Ando, J. F. Beacom, Phys. Rev. Lett., 2005,95: 
171101-1 - 171101-4 

[5] C. W. Akerlof et al., arXiv:astro-ph/0210238vl, 2002, 

[6] T. Totani, A. Panaitescu, The Astrophysical Journal, 
2002,576:120-134 

[7] J. Greiner et al., Astronomy and Astrophysics Supple- 
ment, 1999, 138: 441-442 

[8] Achterberg et al., The Astrophysical Journal, 2007, 
664: 397-410 

[9] M. Kowalski, A. Mohr, Astropart. Phys., 2007, 27: 
533-538 



4 Current Status and Outlook 

The program was approved by the IceCube collaboration at 
the beginning of 201 1 and is running since February, 1 I th . 
Until the midst of May, one alert has been forwarded to 
Swift. The total latency between the neutrino events and 
the first observation by Swift was 90 minutes. All steps 
in the alert chain worked as planned and the event will be 
included in the final analysis of the program. 
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limit on the model is set. The combined limit is 0.22 times the predicted neutrino flux. Finally, the implications for the 
fireball model are discussed. 
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1 Introduction 

Gamma-ray Bursts (GRBs) are prime candidates for the 
production of the highest energy cosmic rays because of 
the enormous energy that is released in such an event [1] 
(O(10 51 — 10 54 crg x 0,/A-k) in gamma rays, where is 
the opening angle of a possible beamed emission). If the 
prime engine accelerates protons and electrons with simi- 
lar efficiencies this would be sufficient energy to account 
for the observed ultra high energy cosmic rays. The ob- 
served gamma-rays would originate from high energy elec- 
tron synchrotron emission and inverse Compton scatter- 
ing, while high energy neutrons would escape the fireball's 
magnetic field and later decay to protons, which would be 
responsible for the high energy cosmic ray flux seen on 
Earth. The observation of high energy gamma-rays con- 
firms the presence of high energy electrons in the fire- 
ball; however, because high energy protons are deflected in 
inter-galactic and the Galactic magnetic fields no direct ob- 
servation of protons from GRBs is possible. Nevertheless, 
if high energy protons are present in the fireball along with 
high energy electrons it is reasonable to assume that pions 
will be produced through pj interactions near the source, 
which would give rise to neutrinos. Guetta et al. [2] gives 
a detailed account of the expected neutrino flux from such 
interactions and is the model that is used for the theoreti- 
cal neutrino prediction in this paper. Previous searches with 
IceCube and other experiments have given null results, with 
the most recent search done in IceCube achieving a 90% 
upper limit that is slightly below the predicted model flux 



[3]. In this contribution, a further improved limit is pre- 
sented, which is then combined with the previous one. 

2 IceCube 

IceCube is a km 3 -scale neutrino detector at the South Pole 
sensitive to TeV-scale neutrinos and above. Construction 
of the detector finished in December, 2010. IceCube de- 
tects Cherenkov light emitted by secondary charged par- 
ticles produced in neutrino nucleon interactions and uses 
that information to reconstruct neutrinos. The finished de- 
tector is made up of 5 160 optical modules (DOMs), with 60 
optical modules placed on each of the 86 strings. The re- 
sults presented here were obtained with the 59-string con- 
figuration of IceCube, which took data from 05/20/09 to 
05/31/10. IceCube is able to detect all known neutrino fla- 
vors; however, in this analysis the focus was on v^. Fur- 
thermore, IceCube is sensitive to the entire sky; however, 
because of the large cosmic-ray muon background in the 
southern sky, this analysis only considers events that were 
reconstructed as coming from the northern sky and con- 
sequently, only GRBs in that part of the sky were anal- 
ysed. In this region, the best sensitivity for can be 
achieved in part because of the good angular resolution for 
muons (0.7° for E v > 10 TeV) and the low background. 
The background consists of mis-reconstructed muons (a 
reducible background) and atmospheric neutrinos (an ir- 
reducible background). Both backgrounds have a softer 
spectrum than the predicted neutrinos from GRBs so event 
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energy information can be used to improve the signal to 
background ratio. 

3 Event Reconstruction 

Events in IceCube are reconstructed by fitting the spa- 
tial and temporal Cherenkov light hit pattern observed by 
the DOMs in a muon event using a maximum likelihood 
method [4] [5] . In the energy range that IceCube is sensitive 
to, neutrinos have sufficiently high energy for the charged 
current interaction between the neutrino and the nucleon 
to be forward and hence the muon and neutrino move in a 
nearly collinear manner, which enables the determination 
of the neutrino direction from the reconstructed muon. The 
shape of the likelihood space near the maximum gives an 
estimate of the reconstruction error of the fit [6]. In ad- 
dition to knowing the direction of the neutrino, knowing 
the energy helps to separate signal from background. The 
stochastic nature of the muon energy loss in the ice, and 
the fact that many tracks originate outside of the detector 
makes it impossible to measure the energy of a muon at 
the neutrino-interaction point directly. Nevertheless, it is 
possible to measure the energy loss rate of a muon as it 
traverses the detector, which is correlated to the energy of 
the muon inside the detector for energies > 1 TeV [7]. 
The energy resolution achieved in this way is 0.3 to 0.4 in 
login (E). 

4 The GRB sample 

During the IC59 data taking period, 105 GRBs were ob- 
served in the northern sky and reported via the GRB Co- 
ordinates Network (GCN) [8]. Of those GRBs 9 had to 
be removed, because IceCube was not taking physics data. 
GRB090422 and GRB090423 happened during 59-string 
test runs before the official start of the IceCube-59 runs and 
were included in the final GRB list as well, which brings 
the final catalog to 98 GRBs. The GRB localization is taken 
from the satellite that has the smallest reported error. The 
start (T s t ar t) and stop (T s t p) times are taken by finding the 
earliest and latest time reported for gamma emission. The 
fluence, and gamma-ray spectral parameters are taken pref- 
erentially from Fermi (GBM), Konus-Wind, Suzaku WAM, 
and Swift in this order. The gamma-ray spectra reported by 
the satellites were used to calculate the neutrino spectra and 
flux as outlined in Appendix A of [2]. The neutrino energy 
spectrum was calculated as a power law with two breaks, 
with the first break corresponding to the break in the pho- 
ton spectrum and with the second break corresponding to 
synchrotron losses of muons and pions (Fig 1). GCN does 
not always report values for all of the parameters used in 
the neutrino spectrum calculation. In that case average val- 
ues are used for the parameters not measured by the satel- 
lites. GRBs are classified into two groups: long soft bursts, 
which are all bursts with a duration longer than 2 seconds 
and short-hard bursts, which are all bursts with a duration 



of less than 2 seconds. Average parameters from [3] were 
used. 




E. [GeV\ 

Figure 1: Neutrino spectra of the GRBs used in this anal- 
ysis. The thin lines represent the individual bursts while 
the solid thick line represents the sum of all bursts. Finally, 
the dashed line shows the Waxman 2003 [9] prediction nor- 
malized to the number of GRBs observed. 



5 Analysis 

The analysis presented here was designed to be sensitive 
to neutrino production from p-f interactions in the prompt 
phase of the fireball. To separate signal from background a 
Boosted Decision Tree [10] was trained. The analysis was 
then optimized for discovery with respect to the Boosted 
Decision Tree score. The optimized value resulted in a fi- 
nal data sample of 85% atmospheric neutrinos and 15% 
miss-reconstructed cosmic ray muons in the off time data 
sample (any events not within ±2 hours of a GRB). An 
unbinned maximum likelihood search [11] was performed 
and each event passing the boosted decision tree cut was 
assigned a probability of being a signal event from a GRB 
or a background event. The final likelihood is the product 
of three PDFs based on the location of an event with re- 
spect to a GRB, the timing information of the event with 
respect to the prompt gamma-ray emission, and the energy. 
The directional signal PDF is a two-dimensional Gaussian: 

PDFf (a?) = e ^+C> (1) 

2tt{(j 1 + <J GRB ) 

where Oi is the directional uncertainty for the i th event 
and cfgrb is the uncertainty of the GRB location as re- 
ported by GCN. \xi — xqrb\ is the angular difference 
between the reconstructed muon direction and the GRB 
location reported by GCN. The background spacial PDF 
was constructed from off-time data, taking into account the 
direction-dependent acceptance of the detector. 

The time PDF is flat over the duration (Tjqq) of the burst 
and falls off smoothly as a Gaussian on either side. The 
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width, a, of the Gaussian is equal to the T\qq of the burst 
with a minimum of 2 s and a maximum of 25 s . 

The third component of the likelihood is an energy PDF. 
In previous analyses, a single energy PDF for the whole 
northern sky was used [3, 11]. However, because the Earth 
is opaque to neutrinos above ^100 TeV, the northern sky 
was split into three zenith regions in order to account for 
this effect. The signal energy PDFs were computed from 
the reconstructed muon energy-loss (dE/dx) from signal 
simulation and averaged over all GRBs in a region. The 
energy background PDF was computed from the dE/dx 
distributions of all off-time data in each region. 

The final likelihood is maximized by varying the assumed 
number of signal events n s and a test statistic A is com- 
puted from the likelihood ratio L(n s = n s )/L(n s = 0), 
where h s is the number of signal events for the maximized 
likelihood. A distribution, A, for the background-only case 
is constructed from off-time data by scrambling it in time 
a sufficient number of times. By comparing the A value 
for the on-time data with the background-only distribu- 
tion a p-value for the measurement is derived, which is a 
measure for the compatibility of the measurement with the 
background-only hypothesis. 

6 Result 

No events were found in the on-time data to be on-source 
(within 10° of a GRB) and on time with a GRB and the 
likelihood maximization yielded A = 0. In total 24 back- 
ground events (not necessarily on source) were expected to 
be in the total time window and 21 were observed (none 
on-source). From the Guetta et al. model [2] 5.8 sig- 
nal events were predicted and a final upper limit of 0.46 
times the predicted flux can be set. This limit includes a 
6% systematic uncertainty. The systematic uncertainty is 
estimated by varying parameters in the signal simulation 
and recomputing the limit, with the dominant factor being 
the efficiency of the DOMs (the uncertainty of the DOM- 
efficiency is <~ 10 %). 
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Figure 2: This plot shows the result of this analysis along 
with the result of the previous analyses. The flux lines from 
the predictions from Guetta et al. [2] and Waxman 2003 [9] 
are shown as well. 



The corresponding model dependent result presented in a 
previous analysis [3] sets a limit of 0.82 of the model flux. 
This limit was obtained using data from the IceCube detec- 
tor in the 40-string configuration. It is possible to combine 
this limit with the limit presented in this paper [3], because 
both analyses obtain a null result. The limits are combined 
by using signal simulation from each analysis and com- 
bining them into one signal simulation data set. From the 
combined signal data set a new limit is calculated by find- 
ing the fraction of total signal flux that would have yielded 
a test statistic that was greater than zero in either analy- 
sis in 90% of the cases. This new fractional signal flux is 
the combined limit and is 0.22 times the flux calculated ac- 
cording to Guetta et al. [2]. Systematic uncertainties were 
handled by combining the worst limit from each analysis 
which makes the combined limit conservative with respect 
to systematic uncertainties. Figure 3 shows the combined 
limit from these two analyses. 
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Figure 3: The combined limit of the IC40+59 analysis is 
shown in addition to the limits and flux predictions dis- 
played in Fig. 2. 



7 Discussion 

Previous results have excluded the neutrino production 
models outlined in [2] and [9] at a level where it may still 
have been explained by statistical fluctuations. This analy- 
sis is able to exclude the models with high confidence and 
if the result is combined with the previous result the model 
in question is strongly disfavored. The caveat is that there 
are parameters in the model for which average values, or 
theoretically calculated values are used, because they are 
not measured (or rarely measured) by the satellites. The 
bulk Lorentz Factor T is one of these values. The lower 
limit on this value is established by pair production argu- 
ments [2], but the upper limit is less clear. Recent pa- 
pers [12, 13, 14] suggest that T can take values of up to 
1000 (316 was used in this analysis as well as in [11, 3]). 
r is an important parameter, because in this analysis the 
GRBs that have the highest neutrino expectation also have 
the highest energy gamma-rays observed by Fermi's LAT 
[15]. Because of pair production arguments [2], this indi- 
cates higher T factors, which implies that the theoretical 
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brightest GRBs in the neutrino sky would be suppressed 
in practice. Another unmeasured parameter that could con- 
tribute to the non-detection of a neutrino flux from GRBs is 
the variability of the observed 7-ray light curve, t V ar- This 
parameter is assumed to be the characteristic time scale be- 
tween the collision of different shock fronts in the GRB 
fireball. Conceptually, if this time is shorter, shock fronts 
will collide more frequently, causing a greater number of 
accelerated particles and therefore more neutrinos. Recent 
limits on t V ar indicate that if t V ar is varied by a factor of 10 
(either higher or lower) UHECR could still be explained as 
originating from GRBs [16]. Therefore, tyar was varied by 
a factor of 10 and the limit was recomputed in incremental 
steps from 0.1 — 10 times the standard t V ar value. In Fig. 
4 the limit of this analysis is plotted as a function of V and 
t V ar- It is also useful to ask, how well IceCube will do in 
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Figure 4: Limit set by the IC59 analysis as a function of 
the tyar multiplication factor and V assuming the Guetta et 
al. model [2]. The gray scale indicates the fraction of the 
model flux that can be excluded at each point of the phase 
space at 90% CL. The thick black line indicates where 1 x 
the model can be excluded at 90 % CL, while the dashed 
lines indicated the standard values used in this analysis. 
The excluded region is the region found to the left and be- 
low the exclusion line. 

the 86-string configuration with respect to constraining this 
parameter space. To get a handle on this the GRB catalog 
from the 59-string configuration was used to estimate the 
86-string sensitivity from Monte Carlo simulations. The 
result for 3-years of 86-string operation is plotted in Fig. 5. 
As seen from the plot, 3-years of IceCube-86 can exclude a 
large portion of the allowed parameter space with the por- 
tion that is not excluded being disfavored by theory. The 
parameters that are treated here are not the only variable 
parameters in the model that are important to the neutrino 
flux. However, the above parameters alone could account 
for the null result seen in IceCube so far. Future observa- 
tions with the completed IceCube detector will be able to 
exclude or confirm GRBs as the major sources of UHECR 
production in a few years. 
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Figure 5: Projected sensitivity of IC86 after 3 years of op- 
eration with respect to T and t V ar- 
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Abstract: The detection of astrophysical point sources of neutrinos is a prime goal of the IceCube neutrino telescope. 
Probable high-energy neutrino sources of interest include transient events such as core-collapse supernovae and gamma 
ray bursts (GRBs). It has been proposed that jets are present not only in supernovae that lead to long GRBs but also 
more frequently in so called choked GRBs that lack a high-energy electromagnetic signature. Choked GRBs may be 
detectable by IceCube's DeepCore subdetector. The transient nature of these events coupled with the angular direction 
and current filtering algorithms should allow strong background rejection. We will present simulations of choked GRB 
signal at trigger level and with preliminary data selection cuts applied. 
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1 Introduction 

Long duration gamma-ray bursts (GRBs) have shown a 
strong association with core collapse supernovae [1]. The 
leading model of GRBs attributes the production of gamma 
rays to Fermi-accelerated electrons in internal shocks of 
relativistic jets driven by the core collapse of the progen- 
itor [2]. These GRBs require heavy progenitors (M > 25 
M Q ) and highly relativistic jets (Lorentz boost factor T > 
100) that break through the surrounding stellar envelope. 
Although long duration GRBs appear correlated with su- 
pernovae, very few (< 10~ 3 ) observed supernovae them- 
selves are associated with GRBs [3]. 

It is conceivable that a large fraction of core collapse 
SNe produce mildly relativistic jets. Unlike GRBs, these 
jets never breach the stellar envelope and are essentially 
'choked' within the progenitor. The inability of the jet to 
break through the envelope could arise from either the en- 
velope itself being more massive than that of a GRB event 
or simply due to a lack of sufficient energy. These choked 
GRBs could be part of a continuum class of astronomical 
objects with long duration GRBs (having highly relativis- 
tic jets) representing the far end of the spectrum. Recently, 
evidence for mildly relativistic jets has been observed in 
supernovae 2007gr [4] and 2009bb [5], as well as in the 
observed asymmetry in the explosions of core collapse su- 
pernovae [6, 7]. This lends credence to the notion that cen- 
tral engines with less relativistic jets might occur more fre- 
quently than observable, fully developed GRBs. 



Both hidden and visible jets can accelerate protons in 
shocks, resulting in the production of neutrinos. Despite 
lacking an elecromagnetic signature like typical GRBs, 
these choked GRB events would still have an associated 
burst of neutrinos that could provide information about hid- 
den jets. One model for this type of event has been pro- 
posed by Razzaque, Meszaros and Waxman [8], and it has 
been extended upon by Ando and Beacom to include kaon 
production [9]. This model will hereafter be referred to as 
RMW/AB. The neutrino spectrum predicted by RMW/AB 
is fairly soft but high in fluence, and should be within reach 
of the IceCube detector (> 100 GeV sensitivity) and the 
DeepCore subdetector (> 10 GeV sensitivity). Due to the 
soft nature of the RMW/AB spectrum, DeepCore will be 
better suited to detecting choked GRBs. 

IceCube is a neutrino detector located at the South Pole 
optimized for neutrino energies on the TeV scale. Finished 
in December of 2010, it detects Cherenkov light emitted by 
secondary charged particles produced in a neutrino nucleon 
interaction. The completed detector is made up of 5160 
optical modules, with 60 optical modules placed on each 
of the 86 strings. These optical modules contain PMTs 
with onboard digitizers and are more succintly referred to 
as DOMs (Digital Optical Modules). The DeepCore subar- 
ray includes 8 densely instrumented infill strings optimized 
for low energies plus 12 adjacent standard IceCube strings. 

Complete PMT waveforms are recorded by the DOMs that 
meet the Hard Local Coincidence (HLC) condition. HLC 
requires hits in a DOM and its nearest or next to nearest 
neighbor in a time window of ± 1 p sec. IceCube also 
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records compact information for Soft Local Coincidence, 
or SLC, hits that do not meet HLC. DeepCore's trigger re- 
quires 3 HLC hits in a time window of 2.5 fj, sec. 

Despite its smaller detector volume, DeepCore's enhanced 
sensitivity to lower energies greatly increases the observ- 
able flux. In addition, the location of DeepCore inside the 
IceCube detector should allow for significant background 
rejection through utilization of IceCube itself as a veto. 
For these reasons this analysis focuses on simulating the 
response of the DeepCore detector under its standard trig- 
gering and filtering, and we calculate the expected event 
count from a sample choked GRB for the fully completed 
DeepCore under the RMW/AB model. 

2 Neutrino Production in Jets 

The RMW/AB model assumes a mildly relativistic baryon- 
rich jet with a bulk Lorentz factor Tb = 3 and an opening 
angle 6j <~ T^ 1 = 0.3. The kinetic energy of the jet is 
set to Ej = 3xl0 51 erg, a typical energy for GRBs. The 
variability timescale of the engine mirrors that of observed 
GRBs as well and is set as t v ~ 0.1 s. Shocks within 
the jet accelerate protons with a spectrum ~ E~ 2 up to 
a maximum proton energy of 2x 10 6 GeV determined by 
the acceleration timescale and radiative cooling. Neutri- 
nos are the product of the kaons and pions produced in p-p 
interactions of the accelerated protons with the stellar enve- 
lope. Energies and densities involved are similar to those in 
neutrino production in the Earth's atmosphere. In the case 
of neutrinos from pion decay, the neutrino flavor flux ratio 
<j> u<! : <j> Ufi : 4> Vt is 0:1:0. Secondary neutrinos from muon 
decays can be ignored here because the muons from pion 
decay are immediately subjected to radiative cooling. As 
for neutrinos from kaon decay, the small flux of v e from 
decay is neglected by RMW/AB. Thus, a flavor flux 
ratio of 0:1:0 is also assumed for neutrinos from kaon de- 
cay. After accounting for vacuum oscillations, the expected 
flavor flux ratio at Earth becomes <~1 :2:2 for both contribu- 
tions. Neutrinos are emitted over a time window of O(^10 
s), set by the star's size (At ~ R*/c). 

The shape of the neutrino spectrum is dependent upon that 
of the mesons from the p-p interactions. Initially, these 
mesons have the same E~ 2 spectrum as the protons, but 
mesons undergo hadronic and radiative cooling before de- 
cay. The result is a meson spectrum with two break en- 
ergies at which the spectrum becomes steeper. The neu- 
trino spectrum will match the meson spectrum, and it can 
be modeled as a doubly broken power law. For a given su- 
pernova at 10 Mpc with Tb = 3, opening angle 6j <~ T^ 1 = 
0.3, and Ej ■ = 3 x 10 51 erg, the spectrum is of the form: 
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F v is the all flavor flux normalization where d<fr v /dE is 
5 x 10~ 2 GeV- 1 cm- 2 (5 x 10~ 5 GeV^cm- 2 ) atE^ for 
pions(kaons). The break energies E^ and E|> 2 ^ denote the 
onset of hadronic and radiative cooling respectively where 
Ei 1} =30 GeV (200 GeV) and E^ 2) =100 GeV (20 TeV) for 
pions(kaons). The neutrino flux from both pion and kaon 
contributions is shown as a function of energy in Fig. 1. 
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Figure 1 : All flavor neutrino flux from pion, kaon and com- 
bined contributions. 



3 Effective Area of DeepCore and Expected 
Events 

We calculate the expected number of observed Ice- 
Cube+DeepCore neutrino events N b s given a flux d<& v /dE 
and detector effective area A e f / by 

N obs = JdEA eff (E)^ (2) 

In order to properly estimate the number of expected events 
from an astrophysical source, the neutrino effective area 
must be calculated through detailed simulation of a bench- 
mark incident flux and the detector hardware. We briefly 
describe the calculation of the effective area. 

The effective area of the detector has been calculated by 
simulating neutrinos in the nearby volume surrounding the 
detector, propagating them, and forcing them to interact 
(preventing the simulation of events that do not interact 
within the volume). These events are re- weighted to reflect 
the probability that the interaction would actually occur. 
All flavors of neutrinos used in this proceeding were simu- 
lated with NUGEN (a modified version of ANIS [10] that 
works with IceCube software). Simulation with NUGEN 
includes several effects including the ice/rock boundary be- 
low the detector, Earth neutrino absorption, neutral current 
regeneration, etc. The flavor flux ratio is taken from the ra- 
tio predicted by the RMW/AB jet model for neutrinos orig- 
inating from both pions and kaons. Vacuum oscillations are 
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included [11], and the ratio of neutrino to anti-neutrino is 
assumed equal for all flavors. The propagation of muons 
within the detector has been simulated with MMC [12]. 
Detection of events is determined by simulating the detec- 
tor response to light produced by the daughter lepton (or 
cascade) of the interacting neutrino. Events are considered 
detected if they activate the standard DeepCore trigger. The 
calculated neutrino effective area under the standard Deep- 
Core trigger is shown for all flavors in Fig. 2. 
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Figure 2: Effective area of the DeepCore detector given 
standard DeepCore SMT3 triggering for all flavor {v e - 
triangle symbols, ^-dashed line, ^ r -dashdot line). The 
effective area has been averaged over the entire sky. 

As Fig. 2 shows, higher energy events are significantly 
more visible than either v e or v T . This can be attributed 
to muons produced outside of the physical DeepCore vol- 
ume that then propagate near or through the ice occupied by 
DeepCore DOMs. At lower energies however, these muon 
tracks become shorter and more closely resemble cascade 
events. 

It should be noted that there are some aspects of the simu- 
lation which are not accurate. Although NUGEN performs 
adequately at higher energies typical of IceCube analyses 
(> 100 GeV), the cross-sections, and consequently interac- 
tion probabilites it predicts, lose accuracy at lower energies 
(between 10-100 GeV). This is particularly true for v T as 
the simulation used does not properly take into account the 
kinematics of the r lepton, and it is likely that the actual 
rate of v T will be appreciably lower. 

Combining the calculated effective areas with the flux pre- 
dicted by the RMW/AB model via Eq. 1 yields an esti- 
mation on N bs for the DeepCore detector. The number 
of expected events by flavor and the predicted background 
rate are listed in Table 1 . 

The result for a reference supernova at 10 Mpc is 
an all-flavor expectation of ^10.5 events in the Ice- 
Cube+DeepCore detector under standard DeepCore trig- 
gering. This event expectation is subject to large variation 
due to uncertainties in the jet parameters of the RMW/AB 
model. After application of the DeepCore filter, which 



Flavor 


Trigger 


Filter 


Preliminary Data Cuts 


Ve 


1.6 


1.5 


1.5 


V V 


4.6 


3.9 


3.3 


V T 


4.3 


3.6 


3.1 


Corsika 


250 Hz 


7 Hz 


1.2 Hz 



Table 1: Event expectation in DeepCore by flavor for 
RMW/AB model choked GRB at lOMpc. Background is 
simulated with CORSIKA, and rates are estimated by tak- 
ing the product of the DeepCore trigger rate and the sim- 
ulated rejection factor at each cut level. Event estimation 
for v T may be overly optimistic due to issues in NUGEN 
simulation at lower energies. 

discards events that show causal relation to hits in the 
non-DeepCore IceCube strings (veto region), the all-flavor 
event expectation is about 9 events. This decrease in events 
is mostly due to the rejection of and v T interactions out- 
side of DeepCore's fiducial volume. Events with the best 
possibility of reconstruction will be those neutrinos that in- 
teract within the DeepCore fiducial volume. Therefore, it 
is of interest to estimate the number of events that actually 
interact inside DeepCore. Examination of the simulated in- 
teraction vertices reveals the number of trigger level events 
originating in DeepCore to be 5.7 (1.1 due to v e , 2.4 due to 
v^, and 2.2 due to v T ). The predicted event spectra for all 
flavors at trigger level is shown as a function of energy in 
Fig. 3. 
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Figure 3: Expected trigger level signal in DeepCore as a 
function of energy for all flavors. 

The plot shows a peak in event expectation at about 40 GeV 
for all flavors. This places most of the expected events be- 
low the typical IceCube threshold and well into the energy 
range of DeepCore. 

4 Atmospheric Muon and Neutrino Back- 
ground 

A major goal of the DeepCore detector is to open up the 
southern sky to analysis to obtain a full Air sr view. In order 
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to do this, extensive steps towards reducing the large atmo- 
spheric muon background must be taken. The expected rate 
of atmospheric neutrinos in DeepCore is about 10 5 events 
per year over 2ir while the atmospheric muon rate in Deep- 
Core is a factor of ^10 6 larger. We are currently develop- 
ing techniques to be used in conjunction with the IceCube 
veto that will allow for a rejection of atmospheric muons 
by a factor > 10 6 while maintaining high signal efficiency. 
For more information on DeepCore and its background re- 
jection capabilities, see reference [13]. 

We have already begun investigating other background re- 
jection methods. Some possible simple cuts include a mod- 
ified version of the current filter, an algorithm using a re- 
construction to search for correlated single detector hits 
(such hits are often cleaned in analysis), and a cut on the ra- 
tio of DOM hits inside and outside of the DeepCore detec- 
tor. Taken in combination, these three additional data cuts 
can reduce background after standard filtering by an addi- 
tional factor of six while maintaining a signal efficiency of 
^94%. The effect of these cuts on the number of expected 
events is shown in Table 1. Any future analyses will re- 
quire a reduction in background to about the atmospheric 
neutrino level (~ 3 mHz). 

5 Discussion 

By simulating a choked GRB in accordance with the 
RMW/AB model, we have predicted the expected neutrino 
event count in the IceCube+DeepCore detector. For a ref- 
erence supernova at 10 Mpc with a bulk Lorentz boost fac- 
tor Tb = 3, opening angle 6j <~ T^ 1 = 0.3, jet energy 
E 3 ■ = 3 x 10 51 erg, and time variability t v ^0.1 s, we ex- 
pect ^10.5 trigger level neutrino events in DeepCore. This 
level of event expectation would make a search for neu- 
trinos in coincidence with known supernova on a distance 
scale ~ 10 Mpc possible. 

One possibile sensitivity enhancement is the expansion of 
DeepCore to a 2-layer IceCube veto (roughly doubling the 
detector volume used). This option would be particularly 
useful in a search for correlated neutrinos from known 
sources. The higher rate of background acceptance brought 
on by expanding the detection volume should be mitigated 
by the increase in signal rate for an overall increase in de- 
tection capability. 

We also expect the event rate predictions to be modified 
when a more accurate simulation of the detector response 
is implemented. Most of this improvement will come from 
the use of superior neutrino cross-section simulation pro- 
vided by GENIE [14]. A version of GENIE compatible 
with IceCube software has recently been developed to sim- 
ulate lower energy neutrinos with much greater accuracy 
than that of the NUGEN based simulation used presently. 
This will greatly improve the simulation of v T events in 
particular, which, due to oscillations, constitute ~ 40% of 
the incident flux. 



One promising method for searching for neutrino bursts 
from choked GRBs is a rolling time window search [15]. 
In this type of search, a time window for bursts is set by the 
characteristic neutrino emission time (At <~ RJc). This 
fixed time window slides across the dataset looking for a 
statistical excess of events. For choked GRBs, the time 
window would be about 10-100 s. Only background events 
falling inside a time window seeing an excess would be 
kept, thus greatly reducing the amount of background ac- 
cepted. An advantage of the rolling search is that it is not 
dependent on any optical observations, allowing it to look 
for photon-dark neutrino sources as expected for choked 
GRBs. 

In addition to searching for choked GRBs, it may also be 
possible for DeepCore to detect neutrinos from high lumi- 
nosity GRBs. A model using parameters inferred from ob- 
servations by Fermi developed by R Meszaros and M.J. 
Rees predicts a neutrino spectrum of luminosity compa- 
rable to the photon component [16]. The model predicts 
a muon neutrino energy spectrum centered around ^12 
GeV. We are curently investigating the event rate expected 
in DeepCore. 

Observations of neutrinos in the DeepCore detector on the 
order of 10-100 GeV in coincidence with supernova would 
be strong evidence for the existence of choked jets from the 
central engine. Such an observation would help to uncover 
the relationship between long duration gamma ray bursts 
and core collapse supernovae, a relationship that is not cur- 
rently fully understood. 

References 

[I] S. E. Woosley and J. S. Bloom, Ann. Rev. Astron. As- 
trophys., 2006, 44: 507-556 

[2] B. Zhang, P. Meszaros, Int. J. Mod. Phys. A, 2004, 19: 
2385. 

[3] E. Berger et al., Astrophys. J., 2003, 599: 408. 

[4] Z. Paragi et al., Nature, 2010, 463: 516. 

[5] A.M. Soderberg et al., Nature, 2010, 463: 513. 

[6] J.R. Maund et al., Astrophys.J (to be published); 

arXiv:0909.2841. 
[7] M. Tanaka et al., Astrophys.J, 2009, 700: 1680 
[8] S. Razzaque, P. Meszaros and E. Waxman, Phys. Rev. 

D, 2003, 68 (083001). 
[9] S. Ando and J.F Beacom, Phys. Rev. Lett., 2005, 95 

(061103). 

[10] A. Gazizov and M. Kowalski, Comput. Phys. Com- 
mun., 2005, 172: 203. 

[II] S. Razzaque, A. Yu. Smirnov, JHER 2010, 3: 31. 
[12] D. Chirkin and W. Rhode (unpublished). arXiv:hep- 

ph/0407075. 

[13] C. Wiebusch et al., in proceedings of the 31st Int. Cos- 
mic Ray Conf., Lodz, Poland, 2009; arXiv:0907.2263. 

[14] C. Andreopoulos et al., Nucl.Instrum.Meth.A, 2010, 
614: 87-104. 

[15] Achterberg A. et al., Astrophys. J., 2007, 664: 397. 
[16] P. Meszaros and M.J. Rees, ApJ, 201 1, 733: L40. 



36 



32nd International Cosmic Ray Conference, Beiiing 201 1 




Neutrino triggered high-energy gamma-ray follow-up with IceCube 

The IceCube Collaboration 1 
1 See special section in these proceedings 



Abstract: We present the status of a program for the generation of online alerts issued by IceCube for gamma-ray 
follow-up observations by Air Shower Cherenkov telescopes (e.g. MAGIC). To overcome the low probability of 
simultaneous observations of flares of objects with gamma-ray and neutrino telescopes a neutrino triggered follow-up 
scheme is developed. This mode of operation aims at increasing the availability of simultaneous multi-messenger data 
which can increase the discovery potential and constrain the phenomenological interpretation of the high energy emission 
of selected source classes (e.g. blazars). This requires a fast and stable online analysis of potential neutrino signals. We 
present the work on a significance based alert scheme for a list of phenomenologically selected sources. To monitor the 
detector and the alert system reliability, monitoring systems have been implemented on different levels. We show data 
from the first weeks of running this system. 

Corresponding authors: Robert Franke 2 (robert.franke@desy.de), Elisa Bernardini 2 (elisa.bernardini@desy.de) 
2 DESY Zeuthen, 15738 Zeuthen, Germany 



Keywords: IceCube, neutrino, NToO, blazars 



1 Introduction 

The major aim of neutrino astrophysics is to contribute to 
the understanding of the origin of high energy cosmic rays. 
A point-like neutrino signal of cosmic origin would be an 
unambiguous signature of hadronic processes, unlike 7- 
rays which can also be created in leptonic processes. The 
detection of cosmic neutrinos is however very challenging 
because of their small interaction cross-section and because 
of a large background of atmospheric neutrinos. Parallel 
measurements using neutrino and electromagnetic obser- 
vations (the so-called "multi-messenger" approach) can in- 
crease the chance to discover the first neutrino signals by 
reducing the trial factor penalty arising from observation 
of multiple sky regions and over different time periods. In 
a longer term perspective, the multi-messenger approach 
also aims at providing a scheme for a phenomenological 
interpretation of the first possible detections. 

The search of occasional flares with a high-energy neutrino 
telescope is motivated by the high variability which char- 
acterizes the electromagnetic emission of many neutrino 
candidate sources. Recent results obtained by the IceCube 
Collaboration [1] indicate that high-energy neutrino tele- 
scopes have reached a sensitivity to neutrino fluxes com- 
parable to the observed high energy gamma-ray fluxes of 
Blazars in the brightest states (e.g. the flares of Markarian 
501 in 1997 [2] and Markarian 421 in 2000/2001 [3]). With 
the assumption that the possibly associated neutrino emis- 



sion would be characterized by a flux enhancement compa- 
rable to what is observed in gamma-rays in such states, neu- 
trino flares could be extracted from the sample of neutrino- 
like events with a reasonable significance. 

These astrophysical neutrinos can be searched for in sev- 
eral ways. Here we present a methods for a neutrino point 
source search that looks for events coming from a restricted 
angular region, which could be identified with a known as- 
trophysical object. Finding neutrino point sources in the 
sky means to locate an excess of events from a particular 
direction over the background of cosmic-ray induced atmo- 
spheric neutrinos and muons. These events might present 
additional features that distinguish them from background, 
for example a different energy spectrum or time structure. 
For sources which manifest large time variations in the 
emitted electromagnetic radiation, the signal-to-noise ra- 
tio can be increased by searching for periods of enhanced 
neutrino emission (a time-dependent search). Of special 
interest is the relation of these periods of enhanced neu- 
trino emission with periods of strong high-energy 7-ray 
emission. However, as Imaging Air Cherenkov Telescopes 
(IACTs) have a small field-of-view and are not continu- 
ously operated such correlation studies are not always pos- 
sible to do after the fact. Therefore it is desirable to ensure 
the availability of simultaneous neutrino and high-energy 
7-ray data for periods of interests. This is achieved by an 
online neutrino flare search that alerts a partner IACT ex- 
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periment when a possible neutrino flare from a monitored 
source is detected. 

Such a Neutrino Triggered Target of Opportunity program 
(NToO) using a list of pre-defined sources was devel- 
oped already in 2006 using the AMANDA array to initiate 
quasi-simultaneous gamma-ray follow-up observations by 
MAGIC [4]. We present here a refined and enhanced im- 
plementation using the IceCube neutrino detector. 
IceCube is a one cubic kilometer neutrino detector oper- 
ating in the glacial ice at the geographical South Pole. It 
consists of 86 strings equipped with 5160 digital optical 
modules (DOMs). Each DOM contains a photomultiplier 
tube to detect Cherenkov light of charged ultra-relativistic 
particles. 



2 Neutrino event selection 

The basis for the neutrino event selection is an on-line fil- 
ter that searches for high-quality muon tracks. The full- 
sky rate of this filter is about 35 Hz for IceCube in its 
2010/2011 configuration with 79 deployed strings. This 
rate is strongly dominated by atmospheric muons. As the 
computing resources at the South Pole are limited one can 
not run more elaborate reconstructions at this rate, so a fur- 
ther event selection has to be done. This so called Online 
Level2 filter selects events that were reconstructed as upgo- 
ing (9 > 80°, 8 — 0° equals vertically down-going tracks) 
with a simple likelihood reconstruction that only takes into 
account the arrival time of the first photon at each Digital 
Optical Module. By requiring a good reconstruction qual- 
ity the background of misreconstructed atmospheric muons 
is reduced. The parameters used to assess the track quality 
are the likelihood of the track reconstruction, the number 
of unscattered photons with a small time residual w.r.t. the 
Cherenkov cone and the distribution of these photons along 
the track. The reduced event rate of approximately 3.6 Hz 
can then be reconstructed with more time intensive recon- 
structions, like a likelihood fit seeded with different tracks 
(iterative fit) and a likelihood-fit that takes into account the 
total number of photo-electrons registered in each module 
(multi-photoelectron fit). Based on this reconstruction the 
final event sample is selected by employing a zenith an- 
gle cut of 9 > 90° for the multi-photoelectron fit and fur- 
ther event quality cuts based on this reconstruction. These 
cuts are optimized to achieve a good sensitivity for flares 
of different time durations. The event selection results in 
a median angular resolution of 0.48 ° for an E~ 2 signal 
neutrino spectrum, the median resolution for events with 
E > 10 6 GeV is < 0.4 °. For each event an angular uncer- 
tainty estimate is calculated. 

The resulting event rate compared to the rate of atmo- 
spheric neutrinos as predicted by Monte Carlo as a function 
of zenith angle can be seen in Figure 1 . 
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Figure 1 : Comparison of the rate of selected events to the 
predicted rate of atmospheric neutrinos based on Monte- 
Carlo for IceCube in its 2010/2011 configuration with 79 
deployed strings. The atmospheric neutrino prediction is 
based on the Bartol conventional flux model and the Nau- 
mov prompt flux model. 

3 The time-clustering algorithm 

The timescale of a neutrino flare is not fixed a-priori and 
thus a simple rolling time window approach is not ade- 
quate to detect flares. The time clustering approach that 
was developed for an unbiased neutrino flare search [6] 
looks for any time frame with a significant deviation of 
the number of detected neutrinos from the expected back- 
ground. The simplest implementation uses a binned ap- 
proach where neutrino candidates within a fixed bin around 
a source are regarded as possible signal events. To exploit 
the information that can be extracted from the estimated re- 
construction error and other event properties like the energy 
an unbinned maximum-likelihood method is under devel- 
opment. 

If a neutrino candidate is detected at time ti around a source 
candidate the expected background N££ is calculated for 
all other neutrino candidates j with tj < ti from that source 
candidate. To calculate N££ the detector efficiency as a 
function of the azimuth angle and the uptime has to be 
taken into account. The probability to observe the multi- 
plet (i, j) by chance is then calculated according to: 



E 



- g bck 



k=N t i 3 -l 

obs 



(1) 



where N obs is the number of detected on-source neutrinos 
between tj and tj. It has to be reduced by 1 to take into 
account the bias introduced by the fact that one only does 
this calculation when a signal candidate is detected. As 
typical flares in high energy gamma-rays have a maximal 
duration of several days we constrain our search for time 
clusters of neutrinos to 21 days. 

If the cluster with the highest significance exceeds a certain 
threshold (e.g. corresponding to 3 a) the detector stability 
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Figure 2: Neutrino flux needed from a given source dec- 
lination to trigger a flare with a significance of 3 a with a 
probability of 50%. The neutrino spectrum is assumed to 
be an unbroken power law with a spectral index of —2. 
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extensive monitoring of the DAQ system and South Pole 
on-line processing. However, most of the information is 
only available with a certain delay after data-taking and 
thus not useful for a follow-up program which requires fast 
alerts. To ensure that alerts are triggered by neutrino multi- 
plets that were detected during stable running conditions a 
simple but powerful stability monitoring scheme has been 
developed. It is based on a continuous measurement of the 
relevant trigger and filter rates in time bins of 10 minutes. 
These rates are inserted into an SQL database at the South 
Pole and are generally accessible a few minutes after the 
respective time bin ended. The rates and ratios of rates 
relevant for the selection of good quality neutrino-induced 
muon tracks are compared to an exponential running 
average of these rates to detect significant deviations. The 
running average is necessary as slow seasonal changes in 
the atmosphere and faster weather changes influence the 
rate of atmospheric muons which dominate the Level-2 
rate. This system was tested off-line on data from IceCube 
in its 59-string configuration and proved to correlate very 
well with the extensive off-line detector monitoring. The 
fraction of data that has to be discarded because it was 
flagged as bad by this method was about 1.6 %. 
To generate a sufficient number of alerts to monitor 
the alert generation and forwarding itself we add 2000 
so-called monitoring sources to the sourcelist (see Section 
5). They are randomly distributed over the northern sky. 
To guarantee blindness for these sky locations the alerts 
for the monitoring sources are generated from blinded 
data events. The blindness is achieved by using the 
previous event time in the transformation from detector 
to sky coordinates for the current event instead of its own 
time. Due to the low event rate on the order of 10~ 3 Hz 
this results in a sufficient random shift of the event right 
ascension. 



Figure 3: Expected number of accidental background alerts 
per year for a source at a declination of 14° as a function of 
the alert threshold expressed in units of standard deviations 
corresponding to a one-sided p-value. 

will be checked and an alert will be sent to a Cherenkov 
telescope to initiate a follow-up observation. Figure 2 
shows the flux needed as a function of declination for a 
neutrino spectrum with a spectral index of —2 to trigger a 
flare with a significance of 3 a with a probability of 50 %. 

To not overwhelm the partner experiment with follow-up 
requests one has to know the number of accidental back- 
ground alerts caused by atmospheric neutrinos. This is 
shown in Figure 3 as a function of the alert threshold. 

4 Stability monitoring 

Data quality is very important for any online alert program 
to minimize the rate of false alerts due to detector or 
data aquisition (DAQ) instabilities. IceCube has a very 



5 Sources 

For a test run of this program we used selection crite- 
ria based on FERMI measurements [5]. For the galactic 
sources we choose sources that were observed in TeV and 
had a FERMI variability index > 15. Blazars were chosen 
according to the following criteria: 

• Redshift < 0.6 

• Fermi variability index > 15 

• Spectral index as observed with FERMI < 2.4 (BL 
Lacs only) 

• FERMI flux 1 - 100 GeV > 1 • 10~ 9 ph cm" 2 s" 1 
(BL Lacs only) 

• FERMI flux 0.1 - 1 GeV > 0.7 ■ 10~ 7 ph cm" 2 s" 1 
(FSRQs only) 
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Figure 4: Angular distribution of events (star symbols) for 
one of the alerts that were generated during a test run of this 
follow-up program from March, 21st 2011 till May, 13th 
2011 (see Section 7). The events contributing to this alert 
were blinded using the procedure described in this paper. 
The weighted average of the contributing events is calcu- 
lated using an event-by-event angular resolution estimator. 
The circle indicates the size of the on-source bin. 



These criteria are motivated by a compilation of different 
hadronic models that provide the guidelines to identify the 
most promising neutrino candidate sources. 22 sources 
(one galactic source, three FSRQs and 18 BL Lac objects) 
were selected according to these criteria in the northern 
hemisphere (5 > 0). 

6 Technical design of the alert system 

After the alerts for this follow-up program are generated at 
the South Pole they are sent to the University of Wiscon- 
sin via the Iridium satellite communication system. This 
low bandwidth connection allows to send short messages 
from the South Pole without any significant delay. Once the 
message arrives in the North it is checked to see whether 
it represents a real alert or a test alert from a monitoring 
source. If it is a real alert, the alert is forwarded to the part- 
ner experiment. Depending on the technical setup this can 
happen e.g. via email or a dedicated socket connection. 

All alerts (real and test) are filled into a database and a 
monitoring web page is updated. Each alert can be re- 
viewed and basic information like the coordinates of the 
contributing events can be inspected. This allows a fast hu- 
man inspection of alerts, even before the full IceCube event 
data arrives in the North. For each generated alert the time 
and space distribution of the contributing events can be in- 
spected (see Figure 4). Furthermore global properties of 
the alerts, like their rate, significance and time length dis- 
tribution are plotted and monitored. 

The total time delay between the time the (latest) neutrino 
event is detected by IceCube and the moment it is for- 
warded to the partner experiment is on the order of several 



minutes (~ lOmin). This time is dominated by the delay 
until the detector rate is available in the database and the 
event processing time in the South Pole system. 

7 Testrun results and Outlook 

The system described here was tested online with the 
IceCube 79-string configuration from March, 21st 2011 
till May, 13th 2011. During the test run neutrino triggers 
were generated online but not forwarded to any IACT 199 
alerts were generated during this test run for all sources 
combined, including the monitoring sources, while 219 
where expected for a 52 day period. Besides statistical 
fluctuations, part of the discrepancy is also due to the 
limited event history available during the first days of 
running the program. 

We plan to run this neutrino triggered high-energy gamma 
follow-up program using IceCube in its final 86-string 
configuration. Several enhancements are possible and 
planned. A maximum-likelihood based significance 
calculation taking into account an event-by-event angular 
reconstruction uncertainty estimation and an energy 
estimation of the event will further improve the sensitivity 
to neutrino flares. 
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Abstract: The observation of a deficit of cosmic rays from the direction of the Moon is an important experimental 
verification of the absolute pointing accuracy of the IceCube detector and the angular resolution of the reconstruction 
methods. This Moon shadow in the downward-going muon flux has been observed with a statistical significance of more 
than 10 sigma in an initial analysis based on a binned counting approach. An unbinned maximum likelihood method 
was developed to reconstruct the shape and the position of this shadow more precisely, to compare the performance of 
different reconstruction algorithms and to verify the correctness of the angular error estimate. 
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1 Introduction 

IceCube [1] is a cubic kilometer scale Cherenkov detec- 
tor at the geographical South Pole, designed to search for 
muons from high energy neutrino interactions. The arrival 
directions of these muons, which are reconstructed with 
0(1°) accuracy, are used to search for point sources of as- 
trophysical neutrinos [2], one of the primary goals of Ice- 
Cube. 

The main component of IceCube is an array of 5160 Dig- 
ital Optical Modules (DOMs) deployed in the glacial ice 
at depths between 1450 m and 2450 m. During construc- 
tion, with the first string of 60 DOMs deployed in January 
2005 and the 86th and final string deployed in December 
2010, the detector already took high quality data. The 
data analyzed and reported here were taken in the 40 and 
59 string configurations, which were in operation between 
April 2008 and June 2010, with a configuration switch in 
May 2009. 

For downward-going directions, the vast majority of the de- 
tected muons do not originate from neutrino interactions, 
but from high energy cosmic ray interactions in the atmo- 
sphere. While these cosmic ray muons are the dominant 
background in the search for astrophysical neutrinos, they 
can be used to study the performance of our detector. In 
particular, we can verify the pointing capability of IceCube 
by studying the shadow of the Moon in cosmic ray muons. 



Cosmic rays at TeV energies propagate through the solar 
system nearly uniformly in all directions. The Moon blocks 
some cosmic rays from reaching the Earth. This creates the 
shadow of the Moon, a relative deficit of cosmic ray muons 
from the direction of the Moon. 

The idea of a Moon shadow was first proposed in 1957 [3], 
and has become an established observation for a number 
of astroparticle physics experiments [4, 5, 6, 7]. Exper- 
iments have used the Moon shadow to calibrate detector 
angular resolution and pointing accuracy [8]. The shift of 
the Moon shadow due to the Earth magnetic field has also 
been observed [9]. 

For an observer at the geographic South Pole, the Moon 
rises and sets once per orbital period of 27.32 days. The 
number of cosmic ray induced muons reaching IceCube 
decreases with increasing declination (i.e. for increasingly 
horizontal directions), since the Earth and the Antarctic ice 
sheet filter low energy muons. Therefore, the shadow of the 
Moon is best observed as far above the horizon as possible, 
i.e. at low declinations. However, the minimum declination 
of the Moon in an orbital period varies slowly over time 
with a period of 18.6 years and is currently increasing. In 
April 2008, 2009 and 2010 the minimum declination of the 
Moon was -26.89°, -25.85° and -24.47°, respectively. 
Fig. 1 shows the energy spectrum of cosmic ray primaries 
that result in one or more muons triggering IceCube. For 
the declination greater than —30°, the energy threshold is 
about 2 TeV. 



41 



D.J. BOERSMA et al. , ICECUBE: SHADOW OF THE MOON 



10- E 



10"' I 



1(F E 




log10(E/GeV) 

Figure 1: The energy spectrum of cosmic ray (CR) pri- 
maries with muons triggering IceCube, as simulated using 
CORSIKA [10]. Dashed: all events; solid: primaries with 

S C r > -30°. 

The Moon shadow in cosmic rays was first observed with 
IceCube using data taken during the first 8 orbital periods 
in the 40-string configuration, using a binned analysis [11]. 
In the analysis using the full data sets from the 40-string 
and 59-string configurations, a log-likelihood based analy- 
sis [12] has now been developed to study the point spread 
function of IceCube for muons. 

The observed Moon shadow can be characterized with the 
following observables: 

• The apparent shift of the Moon shadow from its nom- 
inal position (as computed from the time at which 
each muon event was recorded). A shift of order 0.1° 
is expected due to the Earth's magnetic field. Other 
contributions to a shift could come from e.g. a pos- 
sible bias in track reconstruction or an error in the 
clock used to record the event times. 

• The apparent width and ellipticity of the shadow. 
The apparent radius of the Moon is ~ 0.25°, sig- 
nificantly smaller than the estimated angular resolu- 
tion for muon tracks in IceCube. Hence the width 
of the shadow provides an experimentalal verifica- 
tion of the angular resolution estimate, which could 
for instance be different in zenith and azimuth direc- 
tions. 

• The number of shadowed events should be compat- 
ible with the measured flux of cosmic-ray induced 
muons (at the declination of the Moon) and the solid 
angle subtended by the Moon. Any significantly de- 
viating result would be an indication of a systematic 
error. 



2 Event selection 

The trigger rate from cosmic ray muons was about 1.2- 
1 .3 kHz in the 40-string configuration and close to 2 kHz in 
the 59 string configuration. However, most of those muons 



detected by IceCube travel nearly vertically, and thus they 
cannot have come from directions near the Moon. 

The online event selection is defined as follows: 

• the Moon must be at least 15° above the horizon 

• at least 12 DOMs must register each event 

• at least 3 strings must contain hit DOMs 

• the reconstructed direction must be within 10° of the 
Moon in declination 

• the reconstructed direction must be within 
40°/ cos^) of the Moon in right ascension; 
the cos(<y factor corrects for spherical projection 
effects 

where 5^ denotes the declination of the reconstructed track. 

The online Moon shadow filter was active (i.e., the Moon 
was more than 15° above the horizon) for 7-9 days during 
each 27.3 day orbital period. In that time, between 10M 
and 20M events were selected, depending on the number 
of active installed strings, atmospheric conditions and de- 
tector stability. This is about one percent of all events trig- 
gering IceCube during those days. 

The event sample that passed online selection is subject to 
the same higher level track fitting algorithms as used in the 
searches for point sources of astrophysical neutrinos. The 
track likelihood function used in the fit is based on a sim- 
plified model of the scattering and absorption of light in 
ice [13]. In the offline processing, the track fit is repeated 
using several different seeds. For the majority of the events, 
this leads to a solution which is close to the online fit with a 
slightly improved angular resolution, when studied in sim- 
ulated data. For a fraction of all events, the track fit is am- 
biguous and the iterative fit may yield a completely differ- 
ent direction. 

In the Moon shadow analysis, we characterize each event 
by the zenith angle difference AO = On — 0^ (which is 
equivalent to the the declination difference, thanks to the 
unique geographic location of the detector) and the azimuth 
angle difference A<j> = (<pa — <p^) ■ sin^ between the 
direction of the offline reconstructed track and the nominal 
position of the Moon at the time of the event. 

In the analysis, on-source and off-source subsamples are 
defined using the offline reconstruction. They are again 
defined by an angular window, namely \A0\ < 8° and 
\A(f> + <f> g\ < 8°. Here 4> g = 0° for the on-source sam- 
ple, and (f) ff = ±18° for the left and right half of the off- 
source sample, respectively. The on-source samples for the 
full year data sets of the 40-strings and 59 strings configu- 
rations contain 19M and 22M events, respectively. 

The per-event directional error estimate is derived from the 
variation of the track likelihood function near the solution 
obtained with the track fit [14]. It can be characterized ei- 
ther by the 3 parameters describing the lcr error ellipse, or 
by a single average angular error estimate. In this work, we 
use the latter characterization. 
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The reliability of the directional error estimate was studied 
in simulated data, and simple quality selection criteria were 
developed to ensure that the pull (ratio of real and estimated 
angular error) is on average equal to unity. Moreover, for 
numerical stability, only events with an error estimate in 
the range from 0.075° to 1.5° were accepted. About half of 
the events in the on-source and off-source samples satisfy 
all these criteria. 

3 Likelihood analysis 

An unbinned likelihood analysis was applied to both data 
sets, using an approach similar to the likelihood approach 
taken for the IceCube point source searches [2] . The like- 
lihood for the Moon having shadowed n s events centered 
around x s out of the on-source data sample is expressed as: 

N 

L(x s ,n s ) = ^og{^S i + (l-^)B i ), (1) 

i 

where x s — (AO, A<fi) is the position relative to the nom- 
inal Moon position, n s is the number of signal events, N 
is the total number of events, Si is the signal probability 
density, and Bi is the background probability density. Note 
that Eq. 1 includes no explicit energy-dependent term; this 
a major difference between the IceCube Moon analysis and 
the IceCube point source searches. For the Moon shadow, 
we expect the number of signal events to be negative, as the 
Moon produces a deficit. 

The signal probability density function Si was assumed to 
be Gaussian, with a width given for each event by the esti- 
mated error on the reconstructed position [14]. The back- 
ground probability density function Bi was estimated by 
using the normalized (Moon-centered) declination distri- 
bution obtained from the two off-source regions, and by 
assuming a uniform distribution in (Moon-centered) right 
ascension. 

The likelihood (1) was maximized at every point x s in an 
angular grid around the nominal Moon position, allowing 
the number of "signal" events n s to vary. 

4 Results 

The distribution of the reconstructed number of signal 
events n s is shown in Figures 2 and 3 as a function of 
the offset coordinates of the center of the shadow from the 
nominal Moon position. The shadow of the Moon is ob- 
served as a significant deficit centered at the nominal Moon 
position. 

These results are directly compared with the same distribu- 
tions from the off-source samples, as shown in Fig. 4. The 
distributions of the off-source samples are consistent with 
null shadowing effect from the Moon. The RMS values of 
the n s distributions obtained for the left and right halves 
the off-source are considered as two independent estimates 
of the standard deviation of the background fluctuations. 
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Figure 2: The Moon Shadow from the 40-string configura- 
tion (preliminary). See text for details. 
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Figure 3: The Moon Shadow from the 59-string configura- 
tion (preliminary). See text for details. 
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Figure 4: Fluctuation in n s around the background model 
in one half of the off-source sample for the 59-string data 
set (preliminary). See text for details. 
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40 strings 59 strings agreement between expected and observed shadow depth 

15 14 is a verification of the directional error estimate of the de- 

5734 ± 76 8192 ± 91 fault track reconstruction algorithm as used in the analyses 

5326 ± 544 ± 498 8660 ± 565 ± 681 of the data taken with IceCube in these configurations. 

10-llcr 13-15o- 

0.0° 0.0° 

0.o° 0.0° 5 Conclusions and Outlook 



orbital periods 
expected deficit 
observed deficit 
significance 
6 offset 
4> offset 

Table 1 : Results obtained in the Moon shadow analyses of 
the 40-string and 59-string data sets. Note that the two un- 
certainties given for the observed deficit are the estimates 
for the statistical uncertainty obtained from the left and 
right half of the off-source data samples; the numbers do 
not specify a systematic error. 

The distributions of these values have means compatible 
with zero, as expected for the off-source regions. The RMS 
values of these distributions are now used as two indepen- 
dent estimates of the standard deviation of the background 
fluctuations (see Table 1). 

The number of shadowed events found through likelihood 
maximization are also compared with the expected number 
of shadowed events. The expected number is calculated 
from the average apparent radius of the Moon and the off- 
source flux of downward-going muons from the declination 
of the Moon. 

The depth of the observed shadow is compatible with the 
expected number of shadowed events. The significances 
of the shadows shown in Figures 2 and 3 are calculated to 
be 10-llcr and 13-15ct, respectively. The results for both 
detector configurations are summarized in Table 1 . 

The larger number of strings is the main reason for the 
larger significance of the result with the 59-strings con- 
figuration. With more optical sensors, more events are 
recorded, which are reconstructed with better angular res- 
olution. On the other hand, the minimum declination 
reached by the Moon between April 2008 and May 2009 
was more than a degree less than between May 2009 and 
June 2010. Moreover, the 40-string data sample contains 
livetime from one more orbital period. 

To further confirm our results, the data from the 40-string 
configuration have also been analyzed using a different 
track fit algorithm and a corresponding angular error es- 
timate. Simulation studies of this different reconstruction 
algorithm indicated an average pull of 1.55. Without cor- 
recting for this pull, the Moon shadow analysis resulted in 
a central n s value of —3574 ±434, differing by more than 5 
standard deviations from the expectation: —6373 ± 80. Re- 
doing the analysis with the angular error estimates rescaled 
by a factor of 1 .55 resulted in a fitted n s value compatible 
with expectation. 

The log-likelihood-based analysis relies on the error esti- 
mate of the reconstructed track direction, ranging in this 
data sample from 0.075° to 1.5°. If we were not correctly 
estimating the angular errors, then we would find a shadow 
depth which is significantly smaller than expectation. The 



The shadow of the Moon in cosmic rays has been observed 
with a significance of more than IOct in IceCube data col- 
lected between April 2008 and June 2010. The shadow 
depth is compatible with the expected number of shadowed 
events and has no significant systematic offset. These re- 
sults confirm the pointing capability of IceCube. 

We have started performing an observation of the shadow 
of the Sun. The shadow depth of the Sun should be com- 
parable to that of the Moon, but a larger offset is expected 
due to the solar magnetic field. This offset should be corre- 
lated to the energy of the observed muons. Furthermore, if 
there is a component of high energy antiprotons in cosmic 
rays, then this should result in a faint second shadow with 
the opposite offset. 
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